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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multinlv 

Bv 

'i'o  Obtain 

acres 

4,046.873 

square  metres 

cubic  feet 

0.0283 

cubic  metres 

cubic  yards 

0.76 

cubic  metres 

degrees  (angle) 

0.01745329 

radians 

feet 

0.305 

metres 

inches 

2.54 

centimetres 

metric  tons 

1,000.0 

kilograms 

miles  (US  statute) 

1.609 

kilometres 

ounces  (fluid) 

0.0296 

litres 

pounds  (force) 

4.448222 

newtons 

pounds  (force)- 

second/square  foot 

47.88026 

pascals-second 

square  feet 

0.0929 

square  metres 

yards 

0.9144 

metres 
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SAN  FRANCISCO  CENTRAL  BAY  SUSPENDED  SEDIMENT  MOVEMENT 

SUMMER  CONDITION  DATA  COLLECTION  PROGRAM 
AND  NUMERICAL  MODEL  VERIFICATION 

PART  I :  INTRODUCTION 

Background 

1.  The  dredged  material  disposal  site  near  Alcatraz  Island  is  a  natu¬ 
rally  deep  open-water  site  in  San  Francisco  Bay  (Figure  1) .  This  particular 
site  was  first  designated  in  1894  and  is  the  only  such  remaining  site  in  the 
central  portion  of  San  Francisco  Bay.  This  site  has  been  the  preferred  site 
for  disposal  of  the  sediment  materials  from  the  central  bay  region  for  both 
new  work  and  maintenance  dredging.  The  Alcatraz  disposal  site  was  originally 
designated  because  of  the  depth  and  the  high  local  velocities,  which  were 
assumed  to  disperse  and  carry  the  sediments  out  into  the  Pacific  Ocean. 

2.  Several  dredging  projects  potentially  could  use  the  Alcatraz  dis¬ 
posal  site  in  San  Francisco  Bay.  The  Oakland  Harbor  deepening  project  will 
require  the  removal  of  approximately  7  million  cubic  yards*  of  sediment 
material.  Other  projects  include  the  Richmond  project  (1.5  million  cubic 
yards),  the  J.  F.  Baldwin  III  project  (9  million  cubic  yards),  and  the  Navy 
improvement  dredging  project  (1.5  million  cubic  yards).  In  addition,  annual 
maintenance  dredging  is  expected  to  be  approximately  4  million  cubic  yards. 
The  local  public  and  resource  agencies  are  concerned  that  disposal  of  this 
material,  especially  the  fine-grained  (silt  and  clay)  fraction,  might  cause 
turbidity  problems  that  could  adversely  affect  both  fishing  and  water  quality 
conditions.  According  to  researchers,  phytoplankton  primary  production  is 
light  limited  in  the  bay,  and  excessive  turbidity  could  interfere  within  this 
ecologically  important  process.  Deposition  of  fine-grained  sediments  over 
sandy  substrates  might  affect  benthic  organisms.  In  addition,  low  levels  of 
contaminants  associated  with  some  dredged  material  have  the  potential  to  be 
released  into  the  bay  environment  and  the  corresponding  potential  for  bio¬ 
accumulation  and  cycling  of  any  contaminants . 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  4. 
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Figure  1.  Location  map  of  Central  San  Francisco  Bay 


3.  In  1982,  the  US  Army  Engineer  District,  San  Francisco,  determined 
that  the  dredged  materials  being  disposed  at  the  site  were  not  being  dispersed 
at  the  rate  that  they  were  being  introduced  and  furthermore  that  a  permanent 
mound  was  forming.  In  response  to  this  and  the  choice  of  locating  alternate 
but  more  expensive  disposal  sites,  the  San  Francisco  District  formed  the 
Dredged  Material  Disposal  Management  Program  (DMP) .  The  DMP  was  tasked  with 
both  the  determination  of  the  fate  of  the  dredged  materials  disposed  at  the 
Alcatraz  site  and  the  development  of  numerical  models  to  assess  various  other 
possible  disposal  sites  within  San  Francisco  Bay.  A  request  was  made  that  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES)  provide  assistance  with 
these  tasks .  This  report  is  the  result  of  the  fifth  in  a  series  of  interre¬ 
lated  studies  addressing  the  problem  of  open  bay  dredged  material  disposal. 

4.  There  have  been  four  WES  study  reports  previously  prepared  for  the 

DMP: 


Title 

Report 

Number 

Author 

Date 

"Alcatraz  Disposal  Site 

Investigation" 

Miscellaneous 
Paper  HL-86-1 

Trawle 

and  Johnson 

1986 

"Alcatraz  Disposal  Site  Investiga¬ 
tion:  Report  2,  North  Zone  of 
Oakland  Outer  Harbor  and  Richmond 
Inner  Harbor  Sediments" 

Miscellaneous 
Paper  HL-86-1 

Trawle 

and  Johnson 

1986 

"Alcatraz  Disposal  Site  Investiga¬ 
tion:  Report  3,  San  Francisco 
Bay-Alcatraz  Disposal  Site 
Erodibility" 

Miscellaneous 
Paper  HL-86-1 

Teeter 

1987 

"San  Francisco  Bay:  Modeling  System 
for  Dredged  Material  Disposal  and 
Hydraulic  Transport" 

Technical 

Report. 

HL-88-27 

V .  Pankow 

1988 

The  first  two  reports  dealt  with  the  prediction  of  near  field  behavior  from 
the  instant  of  disposal  through  a  few  hundred  seconds  until  most  of  the  mate¬ 
rial  had  descended  to  the  bed.  A  numerical  disposal  model  was  used  for  this 
purpose.  The  third  study  (Teeter  1987)  performed  laboratory  erosion  tests  on 
deposited  and  remolded  Alcatraz  sediment  beds  to  determine  erodibility.  Re¬ 
sults  of  this  study  were  used  with  field  velocity  data  to  predict  the  capacity 
of  the  site  to  disperse  disposed  materials. 

5.  The  fourth  report  (Pankow  1988)  developed  a  hybrid  model  con'  .sting 
of  a  numerical  component  using  verification  data  and  boundary  conditions  from 
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the  physical  model  of  the  San  Francisco  Bay-Delta  Hydraulic  Model  (SFBM) .  A 
two-dimensional  vertically  averaged  hydrodynamic  model  (RMA-2V)  used  a  compu¬ 
tational  mesh  of  San  Francisco  Bay  in  conjunction  with  the  SFBM  data  to  define 
the  tidal  hydraulics  and  overall  circulation  patterns.  The  Disposal  From 
Instantaneous  Dump  (DIFID)  model  was  then  used  to  simulate  the  short-term  fate 
(1  hr  or  less)  of  materials  disposed  from  barge  or  hopper  dredge  in  a  high- 
resolution  grid  of  the  Alcatraz  dump  site.  Although  the  DIFID  model  can  be 
applied  to  the  short-term  fate  of  materials  disposed  at  the  Alcatraz  site,  it 
cannot  evaluate  the  baywide  movement  of  materials  away  from  the  site. 

Finally,  a  sediment  transport  model  (STUDH)  was  used  to  simulate  the  fate  of 
the  disposed  materials  for  3  to  4  hr  following  discharge  using  the  suspended 
sediment  concentrations  predicted  by  the  DIFID  model  and  the  hydrodynamics 
from  RMA-2V. 


Objectives 

6.  The  major  objectives  of  study  tasks  reported  here  were  to  (a,  gage 
the  seasonal  concentration  and  circulation  of  fine-grained  material  in  Central 
San  Francisco  Bay  (Central  Bay)  during  a  typical  low  freshwater  inflow  period 
of  June  through  October;  (b)  review  and  analyze  selected  previous  studies  on 
currents,  suspended  sediments,  and  material  movement  in  San  Francisco  Bay;  and 
(c)  verify  a  numerical  sediment  transport  model  to  the  circulation  and  concen¬ 
tration  of  natural  suspended  sediments  in  the  Central  Bay  area  during  the  same 
low  freshwater  inflow  condition. 

7.  Subsequent  study  tasks  will  use  the  verified  numerical  sediment 
model  to  predict  the  fate  of  the  dredged  material,  to  provide  information  for 
the  management  of  the  Alcatraz  site,  to  evaluate  environmental  impacts  from 
dredge  disposal  operations,  and  to  relate  and  compare  the  predicted  suspended 
sediment  levels  resulting  from  dredged  material  disposal  operations  to  natural 
levels.  Because  fine-grained  material  historically  constitutes  the  majority 
of  the  dredged  material  disposed  at  the  Alcatraz  site  and  since  the  environ¬ 
mental  concerns  to  be  addressed  concern  turbidity  and  contaminated  dredged 
material,  both  phenomena  being  associated  with  fine-grained  material,  this 
study  emphasized  fine-grained  material  and  not  sand  and  sand  movement  as  bed 
load. 
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Approach 


8.  The  overall  study  approach  was  to  use  a  combination  of  field  data 
collection  and  analysis,  physical  model  results,  and  numerical  model  results 
to  address  meteorologic,  hydrologic,  and  hydrodynamic  conditions  during  the 
typically  dry,  low  freshwater  inflow  periods  of  June  through  October,  and  to 
determine  the  influence  of  these  conditions  on  the  fate  of  fine-grained 
dredged  materials  dispersed  from  the  Alcatraz  site.  The  individual  tasks 
involved  in  the  approach  were  as  follows: 

a.  Plan  and  implement  an  intensive  field  survey  to  obtain  veloc¬ 
ity,  suspended  sediment,  and  related  hydrographic  information. 
The  survey  consisted  of  a  2-week  (spring  to  neap  tide)  period 
of  limited  data  gathering  with  tide  gage  and  automatic  discrete 
water  sampler  operation,  and  a  25-hr  (lunar  day)  intensive 
survey  during  a  near-mean  range  tide. 

b.  Analyze  field  survey  and  other  data  in  the  context  of  suspended 
sediment  behavior,  concentrations,  currents,  and  materials. 

c.  Review  the  published  literature  on  San  Francisco  Bay  currents 
and  sediments . 

d.  Reanalyze  previous  physical  hydraulic  model  test  data  in  regard 
to  Central  Bay  flushing. 

e.  Expand  the  numerical  model  mesh  of  San  Francisco  Bay,  a  devel¬ 
opment  based  upon  previous  WES  numerical  efforts . 

f .  Verify  the  expanded  hydrodynamic  numerical  model  to  SFBM  water 
elevations  and  velocities  and  a  prototype  harmonic  tide. 

g.  Verify  the  numerical  sediment  model  using  the  same  expanded 
mesh  to  the  intensive  field  survey  data  and  limited  historic 
data. 
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PART  II:  DATA  COLLECTION  PROGRAM 


Data  Acquisition 

9.  One  of  the  key  factors  in  successful  numerical  modeling  is  the  qual¬ 
ity  and  quantity  of  the  field  data  used.  Prototype  data  were  used  to  identify 
sedimentation  and  transport  processes  and  geometry  of  the  area,  to  provide 
input  for  the  models,  and  to  provide  a  standard  by  which  model  results  could 
be  compared  (verification) .  Nearly  synoptic  data  were  needed  (all  data  col¬ 
lected  at  the  same  time)  to  develop  a  single,  consistent  set  of  boundary  con¬ 
ditions  to  ascribe  to  all  locations  and  to  establish  relationships  between 
processes  at  the  various  stations.  Long-term  monitoring  data  were  also  needed 
to  assess  processes  at  that  time  scale. 

10.  The  field  data  for  this  study  were  collected  by  WES  personnel  in 
September  1988.  The  prototype  data  acquisition  program  was  performed  over  a 
2-week  period  and  included  an  intensive  25-hr  survey  of  tidal  conditions  dur¬ 
ing  approximately  mean  tide.  During  the  week  preceding  and  the  one  following 
the  intensive  survey,  additional  data,  including  water  levels,  salinities,  and 
suspended  sediment  samples,  were  collected.  During  the  intensive  survey, 
three  to  four  stations  along  four  transects  (ranges)  within  central  San 
Francisco  Bay  were  monitored  over  a  25-hr  period  to  encompass  one  full  tidal 
cycle.  The  intensive  survey  prototype  data  acquisition  program  included  tidal 
elevation,  current  speed  and  direction,  salinity,  and  total  suspended  material 
(TSM)  measurements. 


Station  Locations 


11.  The  program  was  designed  to  obtain  an  overview  of  the  suspended 
sediment  and  salinity  conditions  in  Central  Bay  during  various  wind  and  tidal 
conditions.  A  total  of  5  tide  gage  and  15  point  sample  station  locations 
(Figure  2)  were  established  for  sampling  prior  to,  during,  and  after  the 
intensive  survey,  with  adjustments  made  in  the  field  to  suit  actual  condi¬ 
tions.  For  the  intensive  survey,  four  transects  with  three  or  four  stations 
each  were  established  (Figure  3).  The  final  station  and  equipment  locations, 
along  with  the  descriptions  of  the  equipment  used  and  sample  data  sheets,  will 
be  presented  in  subsequent  sections . 
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Equipment  and  Collection  Methods 


12.  This  section  of  the  report  briefly  describes  the  equipment  used  in 
gathering  the  data.  A  more  specific  description  of  deployment  and  location, 
along  with  other  pertinent  data,  will  be  given  later  in  this  report. 

Tidal  elevations 

13.  Tidal  elevations  were  measured  with  Fischer  and  Porter  1550 
punched-tape  water  level  recorders  (Figure  4).  The  system  consisted  of  a 
stilling-well-contained  float  connected  to  a  recording  device  by  a  wire  rope. 
Elevations  were  recorded  to  the  nearest  0.005  ft.  The  float  elevation  was 
punched  on  16-channel,  foil-backed  paper  tape  every  15  min  by  a  timer- 
activated  recording  mechanism.  The  4-in.-diam  stilling  well  was  designed  to 
attenuate  short-period  fluctuations  in  surface  elevation  of  less  than  1  min. 
Bench  tests  of  the  timers  have  shown  them  generally  accurate  to  ±2  min  per 
month  except  for  infrequent  malfunctions  that  can  cause  larger  time  errors. 
Relative  water  level  accuracy  is  affected  by  the  temperature  of  the  water, 
float,  and  supporting  wire  and  salinity  changes  of  the  water  inside  the  well. 
Relative  accuracy  is  considered  to  be  within  0.1  ft  (Coleman  et  al.  1988). 
Over-the-side  equipment 

14.  Three  of  the  four  survey  boats  were  outfitted  with  over-the-side 
equipment  consisting  of  a  current  meter,  direction  indicator,  wire  rope  and 
suspended  weight,  remote  readout  device,  and  a  support  frame  (Figure  5).  Cur¬ 
rent  velocities  were  measured  using  a  Gurley  Model  665  velocity  meter  with  a 
vertical  axis  cup-type  impeller  in  conjunction  with  a  Magnesyn  directional 
indicator.  This  assembly  was  equipped  with  a  streamlined  weight  that  held  the 
sensors  in  a  vertical  attitude  facing  into  the  flow.  A  winch  and  1/8-in.  wire 
rope  supported  by  a  portable  frame  was  used  to  raise  and  lower  the  assembly 
(Figure  5).  An  indicator  on  the  winch  showed  the  sensor  depth  below  the  water 
surface.  The  Gurley  meters  have  threshold  speeds  of  less  than  0.2  fps  and  an 
accuracy  of  ±0.1  fps  for  speeds  less  than  1  fps  (±5  percent  for  speeds  greater 
than  1  fps) . 

15.  Water  samples  for  salinities  and  total  suspended  solids  were  taken 
with  a  12-v  d-c  pump  using  50  ft  of  1/4-in.  ID  plastic  tubing  attached  to  the 
current  meter  support.  Pumps  and  tubing  were  flushed  with  three  system  vol¬ 
umes  before  individual  samples  were  drawn.  The  samples  and  current  measure¬ 
ments  were  taken  at  3  ft  below  the  surface  depth,  middepth,  and  2  ft  above  the 
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bottom  at  depths  less  than  35  ft,  and  at  five  increments  at  greater  depths. 
Pumped  samples  were  stored  in  8-oz  plastic  bottles  and  kept  cool  (4°  C)  for 
shipment  back  to  WES. 

Additional  equipment  used 

16 .  The  fourth  survey  boat  was  equipped  with  an  InterOceans  Model  S-4 
velocity  meter  and  2-H  water  bottle  samplers.  This  boat  was  used  on  the 
Golden  Gate  range  (Range  2  in  Figure  3)  to  obtain  data  at  depths  of  up  to 

200  ft.  As  the  velocity  meter  was  lowered  in  the  water,  water  bottle  samplers 
were  attached  to  the  hydrographic  wire.  When  the  velocity  meter  was  at  the 
greatest  depth,  a  messenger  weight  was  sent  down  the  wire  to  collect  a  series 
of  water  samples. 

17.  The  InterOceans  Model  S-4  velocity  meter  (Figure  6)  is  another  type 


Figure  6.  InterOceans  S-4  velocity  meter 

of  over-the-side  instrument  that  measures  both  current  velocity  and  direction. 
The  10-in. -diam  sphere  was  used  in  a  tethered  configuration  so  that  the  vari¬ 
able  depth  data  could  be  gathered.  The  device  measured  the  current  using  an 
electromagnetic  field  to  sense  current  induced  by  the  movement  of  water 
through  the  field.  A  microprocessor  coupled  with  a  compass  computed  refer¬ 
enced  current  vectors,  which  were  then  recorded  into  the  solid-state  memory. 
The  data  in  the  memory  were  periodically  loaded  onto  3-1/2-in.  diskettes  using 
a  lap-top  computer  (Figure  6),  and  returned  to  WES  for  Final  data  reduction. 
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Synopsis  of  the  Collection  Program 


18.  Installation  of  sampling  equipment  was  initiated  on  30  August  1988. 
As  described  in  paragraph  10,  during  the  week  preceeding  and  the  week  follow¬ 
ing  the  intensive  survey,  point  samples  were  collected  at  the  15  stations 
shown  in  Figure  2  at  about  2-day  intervals.  A  total  of  six  sample  sets  were 
collected,  at  arbitrary  times  usually  near  slack  water.  Tide  gages  were  ser¬ 
viced  over  the  same  2-week  period. 

19.  The  intensive  survey  used  four  boats  to  collect  data  from  stations 
along  four  separate  ranges.  Arrangements  were  made  with  the  San  Francisco 
District  Operations  Branch  to  use  the  45-ft-long  vessel,  GRIZZLY,  in  the 
larger  open  water  areas  of  the  bay.  Each  range  included  three  stations  with 
the  exception  of  Range  4,  which  had  four  stations  (Figure  3).  Sampling  lasted 
about  25  hr  to  include  one  full  tidal  cycle.  Each  station  was  sampled  approx¬ 
imately  once  every  hour.  This  survey  design  provided  synoptic  data,  which, 
when  compiled,  provided  the  means  for  verification  of  a  numerical  model.  The 
intensive  survey  began  at  0800  PST  on  7  September  and  the  final  sampling  was 
completed  at  1000  PST  on  8  September. 

Actual  Conditions  During  the  Survey 

20.  The  weather  conditions  during  the  overall  survey  period  were  cool, 
breezy,  and  partly  overcast.  The  worst  sea  conditions  were  experienced  at  the 
stations  adjacent  to  the  Golden  Gate  Bridge.  In  that  area,  late  afternoon 
winds,  waves,  and  current  speed  and  commercial  shipping  traffic  made  maneuver¬ 
ing  between  stations  and  station  keeping  difficult.  Sampling  boat  drift  was 
rapid  and  erratic.  For  this  reason,  reliable  current  measurements  could  not 
be  made  at  Range  2.  At  Range  1  near  Richmond  Bridge,  an  anchor  line  fouled 
the  sampling  boat's  propeller.  Despite  the  efforts  of  divers,  the  boat's 
maneuverability  could  not  be  restored,  and  only  one  station  (1C)  was  sampled 
over  the  entire  survey. 

Laboratory  Analysis  and  Data  Reduction 

21.  Laboratory  testing  included  salinity,  TSM,  and  fall  velocity  mea¬ 
surement  of  suspended  sediment. 
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a.  Salinities  were  measured  in  the  laboratory  with  an  AGE  Instru¬ 
ment's  model  2100  salinometer.  Accuracy  of  the  instrument  is 
±0.003  ppt.  Standard  seawater,  known  to  better  than 

±0.001  ppt,  was  used  to  calibrate  the  instrument. 

b.  TSM  was  determined  by  filtration  and  gravimetric  analysis  of 
the  samples.  Nuclepore  polycarbonate  filters  with  0.40-/X  pore 
size  were  used.  A  vacuum  system  (8-lb  vacuum  maximum)  was  used 
to  draw  the  samples  through  preweighed  filters.  After  the 
filters  and  holders  were  washed  with  distilled  water,  the  fil¬ 
ters  were  dried  at  105°  C  for  1  hr  and  reweighed.  Finally,  the 
TSM  were  calculated  based  on  the  weight  of  the  filter  and  the 
volume  of  the  filtered  sample. 

c.  Fall  velocities  were  analyzed  on  a  series  of  composites  made  on 
samples  from  middepth  automatic  samplers  from  near  the  Richmond 
Bridge.  Tests  were  performed  in  a  10-cm-diam  by  2-m-high  clear 
tube  using  the  pipette  method. 

22.  Data  reduction  tasks  included  both  transcribing  and  reformatting 
raw  data,  and  entering  data  obtained  from  both  laboratory  testing  and  over- 
the-side  equipment. 


Explanation  and  Selected  Samples  of  the  Data 


23.  The  following  sections  present  and  describe  selected  data  taken 
during  the  survey.  In  general,  most  of  the  data  gathered  were  acceptable; 
however,  as  in  any  large  data  collection  effort,  equipment  malfunction 
occurred.  Additional  stations  were  included  initially  to  compensate  for  a 
possible  malfunction  of  instrumentation.  Due  to  the  volume,  it  was  impracti¬ 
cal  to  include  all  of  the  data  collected  in  this  report;  therefore,  only  a 
representative  sample  is  presented.  The  data  are  on  computer  file  at  WES  for 
future  reference. 

Tide  data 

24.  Figure  7  shows  example  tide  data  tabulation  and  Figure  8  shows  a 
plot  of  the  same  data.  In  the  tabulation,  a  banner  is  followed  by  a  data 
listing.  Each  line  lists  month,  day,  year,  and  sample  time  for  the  first  data 
point.  The  fifth  field  starts  the  actual  elevation  data,  in  hundredths  of 
feet.  Vertical  datums  were  not  established  for  the  tide  gages. 

Over-the-side  data 

25.  Figures  9  through  11  show  plots  of  current  speed,  salinity,  and 
TSM,  respectively.  Each  shows  surface,  middepth,  and  bottom  time-series  data. 
Figure  12  shows  the  depth-  and  tidal-averaged  TSM  values  for  the  intensive 
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Figure  7.  Example  tide  data  tabulation 
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Figure  8.  Example  tide  data  plot 

survey.  Figure  13  shows  the  depth-  and  tidal-averaged  salinity  values  for  the 
intensive  survey. 

Point  TSM  data 

26.  Figure  14  shows  averages  and  standard  deviations  for  six  sampling 
periods  made  between  28  August  and  23  September  1988. 
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Figure  13.  Depth-  and  tidal-averaged  salinities  for  the  intensive  survey 
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PART  III:  INITIAL  PROCESS  ASSESSMENT 


27.  The  long-term  fate  of  fine-grained  materials  dispersed  from  the 
Alcatraz  disposal  site  depends  on  sediment  interactions  with  flows  and  indi¬ 
vidual  particle  sojourns  through  the  Central  Bay  area.  This  part  summarizes 
initial  interpretations  of  field  data  related  to  dry  weather  seasonal  long¬ 
term  sediment  transport  in  Central  Bay.  The  field  data  described  in  the  pre¬ 
vious  part  will  be  further  analyzed  here.  Some  previous  studies  that  have 
quantified  some  relevant  aspects  of  transport  processes  will  also  be 
summarized. 

28.  Studies  on  the  movements  of  natural  sediments  can  be  useful  to 
assess  the  effects  of  dredged  material  dispersion  from  the  Alcatraz  disposal 
site.  Measures  or  parameters  are  sought  that  will  help  answer  the  following 
questions : 

a.  What  is  the  local  flushing  rate  for  the  sediments  dispersed 
from  the  Alcatraz  site  in  Central  Bay?  (How  fast  do  sediments 
move  away  from  the  site,  and  therefore,  what  is  the  local 
buildup  of  suspended  sediments?) 

b.  Where  are  dispersed  sediments  transported?  (What  is  the  dis¬ 
tribution  probability  of  the  dispersed  materials?) 

c.  What  is  the  fate  of  the  dispersed  sediments?  (How  much 
deposits  in  or  is  transported  through  various  segments  of  the 
Central  Bay  area?) 

29.  The  present  study  estimated  field  sediment  fluxes  entering  and 
leaving  Central  Bay  by  direct  observation  of  suspended  transport  over  a  tidal 
cycle.  The  sampled  net  transport  was  generally  in  the  landward  direction,  but 
net  transport  was  a  small  fraction  of  the  tidal  transport.  It  will  be  shown 
in  this  part  that  the  dominant  transport  mechanism  was  tidal  pumping,  and  that 
near-bed  currents  are  very  important  in  this  process.  Previous  studies  have 
examined  tidal-averaged  current  and  bed  form  patterns.  Those  studies  have 
shown  variability  in  Central  Bay,  with  some  areas  being  flood  dominated  with 
respect  to  maximum  near-bed  currents,  and  other  areas  being  ebb  dominated. 

30.  Sediment  budget  approaches  based  on  changes  in  bay  volumes  or  input 
and  output  accounting  have  also  been  used  to  assess  long-term  sedimentation  in 
San  Francisco  Bay.  These  studies  have  been  concerned  with  annual  or  longer 
time  scales.  Since  the  purpose  here  is  to  examine  only  dry  weather  seasonal 
long-term  transport,  annual  sediment  budgets  will  not  be  discussed. 

31.  The  reader  is  referred  to  the  bibliography  in  Appendix  A  for 
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specific  articles  on  San  Francisco  Bay  sediments.  This  bibliography  includes 
references  on  baywide  sedimentation,  including  such  topics  as  the  influence  of 
hydraulic  mining  in  the  1800 's  on  riverborne  sediment  to  the  bay,  the  use  of 
various  sediment  tracers,  and  global  sediment  budgets  for  the  bay.  However, 
this  bibliography  is  not  intended  to  be  exhaustive. 

Bavwide  Suspended  Sediment  Behavior 

32.  In  the  San  Francisco  Bay  system,  sedimentation  and  suspended  sedi¬ 
ments  have  received  only  limited  attention  compared  to  other  parameters  and 
hydrographic  variables  such  as  water  levels,  currents,  and  salinity.  Most 
field  surveys  have  either  not  included  sediment  and  suspended  sediment  mea¬ 
surements  or  have  not  emphasized  such  measurements.  The  complexity  of  sus¬ 
pended  sediment  movement  in  an  estuarine  environment  makes  quantitative  study 
of  sediment  movement  a  difficult  and  costly  undertaking.  However,  ample  evi¬ 
dence  does  exist  to  provide  qualitative  understanding  of  sediment  and  sus¬ 
pended  sediment  behavior. 

33.  General  agreement  exists  among  most  researchers  regarding  the 
qualitative  behavior  and  characteristics  of  riverine  sediment  loading  to  the 
estuary  including  aggregation,  deposition,  suspension,  erosion,  and  circula¬ 
tion  of  sediments  in  the  San  Francisco  Bay  system.  To  provide  background  to 
baywide  suspended  sediment  behavior,  a  synopsis  of  some  of  the  available 
literature  on  general  baywide  sediment  is  provided. 

34.  In  a  study  of  sedimentation  in  Mare  Island  Strait,  an  artificially 
deepened  portion  of  the  lower  Napa  River  at  Carquinez  Strait,  Einstein  and 
Krone  (1961)  reported  on  the  seasonal  sediment  pattern  in  the  San  Pablo  Bay- 
Carquinez  Strait  study  area.  Einstein  and  Krone,  through  field  observations, 
determined  that  despite  higher  suspended  sediment  concentrations  in  bay  waters 
in  the  winter,  much  of  the  deposition  in  Mare  Island  Strait  occurred  during 
the  summer.  They  concluded  that  deposition  during  the  winter  occurs  in  the 
shallow  bays  and  mud  flats.  The  almost  daily  summer  patterns  of  high  mid¬ 
morning  to  late  afternoon  winds  produce  wave  action  that  resuspends  these 
winter-deposited  sediments.  Once  the  sediments  are  resuspended,  ambient  cur¬ 
rents  move  them  to  relatively  quiescent  areas  such  as  Mare  Island  Strait, 
resulting  in  deposition  in  these  areas. 

35.  These  early  observations  by  Krone  and  Einstein  have  been 
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substantiated  by  subsequent  research  and  field  data  evaluation,  and  a  baywide 
sediment  transport  pattern  similar  to  the  Mare  Island  Strait  results  emerges 
from  the  literature.  See,  for  example,  Conomos  and  Peterson  (1976),  Krone 
(1976,  1979),  and  US  Army  Engineer  District  (USAED) ,  San  Francisco  (1977). 

36.  Approximately  80  to  85  percent  of  the  riverborne  sediment  to 
San  Francisco  Bay  originates  from  the  Sacramento-San  Joaquin  River  basin  and 
over  80  percent  of  this  riverborne  sediment  is  contributed  during  the  period 
of  high  freshwater  inflow  in  winter  and  early  spring.  Furthermore,  the  mea¬ 
sured  suspended  sediment  concentrations  reported  in  Conomos  and  Peterson 
(1976)  indicate  a  strong  pattern  of  maximum  suspended  sediment  levels  in 
San  Pablo  and  Suisun  Bay  areas  (the  null  zone  is  normally  located  in  either  of 
these  bays),  the  lowest  levels  near  the  Golden  Gate  in  Central  Bay,  and  some¬ 
what  elevated  levels  in  South  San  Franciso  Bay  (South  Bay) .  The  seasonal  pat¬ 
terns  reported  by  all  researchers  was  riverborne  sediment  supply  and  deposi¬ 
tion  during  the  winter  and  sediment  resuspension,  transport,  and  redeposition 
during  the  summer. 


Tidal  and  Net  Fluxes  of  TSM 


37.  To  evaluate  the  transport  at  the  boundaries  of  Central  Bay,  tidal 
and  net  TSM  fluxes  were  computed  using  the  intensive  survey  data  from 
Ranges  1-3.  Not  all  data  sets  contained  sufficient  information  with  which  to 
make  the  calculation  using  the  most  accurate  method.  For  this  reason,  com¬ 
putations  were  made  using  several  methods  for  comparison.  Tidal  suspended 
sediment  fluxes  were  calculated  by  three  methods  (direct,  decomposition,  and 
tidal  prism).  Results  are  presented  in  Table  1.  The  most  accurate  method  is 
the  direct  method.  The  direct  method  calculated  average  instantaneous  fluxes 
as  the  product  of  currents,  TSM  concentrations,  and  cross-sectional  areas  cor¬ 
rected  for  instantaneous  tide  height.  These  results  were  then  integrated  over 
the  lunar  sampling  day.  This  method  calculated  vertical  averages  by  weighting 
the  middepth  values  twice  that  of  the  surface  and  bottom  values.  The  values 
used  for  the  surface  and  bottom  were  actually  measured  2  ft  down  from  the  sur¬ 
face  and  2  ft  up  from  the  bottom,  respectively.  Only  Range  3  had  a  complete 
data  set.  Current  and  TSM  data  from  sta  1C  were  applied  to  the  entire  cross 
section  for  Range  1,  because  only  partial  data  were  collected  at  sta  1A  and 
IB.  Average  TSM  data  for  the  beginning  of  the  intensive  survey  were 
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Table  1 

TSM  Conditions  and  Fluxes  (Positive  °  Landward') 


TSM  Flux,  1,000  metric  tons/lunar  day 
_ for  Method* _ 


Ranee 

Tidal 

Phase 

TSM 

me/f 

Direct 

Deconmosition 

Tidal 

Prism 

1 

Ebb 

44.9 

-48 

_ 

-58 

(Richmond 

Flood 

37.5 

39 

- 

49 

Bridge) 

Net 

-9 

-37 

-9 

2 

Ebb 

14.9 

-** 

-** 

-41 

(Golden 

Flood 

16.3 

- 

— 

45 

Gate) 

Net 

- 

- 

4 

3 

Ebb 

33.7 

-18 

-29 

(Oakland 

Flood 

36.1 

19 

- 

32 

Bay 

Net 

1 

-4 

3 

Bridge) 


*  See  paragraphs  37-40  for  explanation. 

**  Unreliability  of  current  measurements  at  sta  2  prevents  evaluation. 

comparable  at  sta  1A,  IB,  and  1C.  However,  the  validity  of  the  flux  calcula¬ 
tion  at  Range  1  is  questionable. 

38.  The  decomposition  flux  calculation  method  averaged  instantaneous 
flux  per  unit  areas  over  the  lunar  day.  Various  flux  components  were  calcu¬ 
lated  as  described  later  in  paragraphs  41-43  and  summed  to  a  total.  The  total 
average  flux  per  unit  area  multiplied  by  the  average  cross  section  and  inte¬ 
grated  over  the  lunar  day  provided  the  decomposition  method  estimate  of  net 
TSM  flux  (Table  1) .  Ebb  and  flood  fluxes  are  not  available  from  this  method. 

39.  A  tidal  prism  method  was  also  used  to  calculate  TSM  fluxes  from  ebb 
and  flood  tidal  prisms  and  average  TSM  concentrations .  This  was  the  only  pos¬ 
sible  method  for  Range  2,  where  current  data  were  lacking.  Surface  areas  for 
the  Central  and  South  Bays  were  estimated  by  planimetering  National  Ocean  Ser¬ 
vice  (NOS)  nautical  charts.  Total  surface  area  and  tidal  prism  estimated  by  a 
cubature  method  for  the  San  Francisco  Bay  system  were  calculated  for  an  aver¬ 
age  tide  (NOS  1985).  The  total  tidal  prism  for  the  bay  system  for  the  survey 
was  estimated  by  tide-range  ratio.  Then,  tidal  prisms  passing  through  the 
various  survey  ranges  were  estimated  for  the  observed  tide  by  subtracting 
various  tidal  prism  segments  from  the  total.  Areas,  widths,  and  tidal  volumes 
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passing  survey  ranges  during  the  intensive  survey  are  given  in  the  following 
tabulation: 


Ranee 

Gross-Sectional 
Area .  m2 

Width,  m 

1 

53,800 

5,940 

2 

117,000 

3,750 

3 

66,100 

4,500 

Tidal 

Volumes . 

109 

10  cu  m 

Flood 

Ebb 

Average 

1.30 

1.28 

1.29 

2.75 

2.71 

2.73 

0.88 

0.87 

0.87 

40.  There  was  a  substantial  variation  between  direct  and  tidal  prism 
estimates  of  tidal  fluxes,  possibly  caused  by  an  overestimation  of  tidal 
prisms.  Tidal  prism  flux  results  averaged  44  percent  greater  than  those  of 
the  direct  method.  Range  3  had  the  best  data  coverage  for  any  range.  The 
direct  method  estimate  of  1,000  metric  tons  landward  net  flux  is  only  a  small 
percentage  of  the  18,000-19,000  metric  ton  tidal  flux.  The  tidal  prism  method 
gave  similar,  though  larger,  values  for  flux  at  Ranges  1  and  3.  Thus,  tidal 
flux  at  Range  2  adjusted  for  the  tidal  prism  method's  overestimation  was  prob¬ 
ably  about  30,000  metric  tons.  The  decomposition  method  result  at  Range  3  was 
not  consistent  with  other  methods,  but  differences  with  the  direct  method  were 
still  small  in  the  absolute  sense. 

Flux  Components  for  TSM 

41.  A  flux  decomposition  method  was  used  to  examine  transport  compo¬ 
nents  at  individual  stations,  as  well  as  to  estimate  total  tidal  flux  as  de¬ 
scribed  earlier.  Tidal  fluxes  of  TSM  were  decomposed  based  on  the  following 
equation: 

Flux  of  C  =  A(UoCo  +  UiCi  +  UovCov  +  UivCiv)  (1) 

where 

C  *=  concentration 
A  -  cross-sectional  area 
overbar  -  depth  average 
U  =  velocity 

o  =»  steady  (tidal-averaged) 
i  -  fluctuating  component 

v  **  vertical  deviations  from  the  depth  mean 
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Tidal  average  values  for  each  of  the  four  terms  inside  the  parentheses  were 
calculated  from  intensive  survey  data.  The  first  term  is  mean-flow  transport. 
The  second  and  fourth  terms  are  tidal  pumping  components .  The  third  term  is 
transport  by  vertical  circulation. 

42.  TSM  flux  decomposition  results  are  shown  in  Table  2.  Also  included 
in  Table  2  are  the  principal  salinity  flux  component,  UoSo  ,  and  salinity  and 
TSM  flux  totals.  Note  that  all  other  salinity  flux  terms  were  of  minor  impor¬ 
tance,  and  that,  except  for  mean-flow  transport,  salinity  conditions  were 
nearly  balanced.  Analysis  of  several  weeks  of  salinity  samples  also  indicated 
that  salinity  conditions  were  steady.  The  mean  flow  component,  Uo  ,  does  not 
always  represent  a  true  tidal  residual  flow  and  can  be  influenced  by  sampling 
bias . 

43.  The  most  important  or  dominant  flux  component  was  UivCiv  ,  the 
tidal  pumping  component  involving  vertical  deviations  in  instantaneous  or 
fluctuating  fluxes.  Near-bed  TSM  fluctuations  were  the  greatest,  as  shown  in 
the  example  (Figure  11) .  The  pattern  of  net  TSM  flux  suggested  by  the  decom¬ 
position  totals  is  indicated  in  Figure  15. 


Table  2 

Station  Average  Flux  Components  (Positive  -  Upstream) 


Salinity  Flux 

fps  x  ppt*  _ TSM  Flux,  fps  X  mg / SL 


Range 

Station 

UoSo 

Total 

UoCo 

UiCi 

UovCov 

UivCiv 

Total 

1 

C 

0.00 

1.07 

0.0 

-9.13 

-0.06 

-28.28 

-37.47 

3 

A 

4.23 

4.68 

4.79 

0.77 

0.21 

-8.07 

-2.30 

B 

6.37 

6.36 

7.68 

2.41 

1.73 

-28.30 

-16.48 

C 

10.90 

10.80 

8.68 

1.59 

-0.13 

0.60 

10.74 

4** 

A 

-0.35 

-0.59 

-0.25 

7.71 

-0.13 

-33.16 

-25.84 

B 

3.62 

3.47 

2.36 

1.68 

0.15 

-15.92 

-11.72 

C 

3.34 

3.55 

2.39 

1.80 

0.05 

-3.36 

0.88 

D 

4.54 

4.67 

3.88 

2.34 

0.09 

1.26 

7.57 

*  1  fps  =  0.305  m/sec. 

**  Positive  =  south. 
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SAN  FRANCISCO 


122° 


122°-20'W 


Figure  15.  Relative  TSM  net  fluxes  according  to  station  flux 

decomposition  data 


Trace  Metals  on  TSM 


44.  Certain  trace  metals  were  analyzed  on  the  suspended  material  col¬ 
lected  on  filters  during  TSM  sample  analysis.  All  filters  from  each  range 
were  divided  into  ebb  and  flood  composites  for  analysis.  Thus  each  composite 
was  made  up  of  roughly  100  filter  subsamples.  Rax;  trace  metals  on  TSM  and 
normalized  trace  metals  on  TSM  are  given  in  Tables  3  and  4,  respectively.  The 
trace  metals  manganese  (Mn) ,  zinc  (Zn) ,  copper  (Cu) ,  and  nickel  (Ni)  were 
normalized  with  reference  to  iron  (Fe). 

45.  Trace  metals  on  suspended  matter  can  indicate  their  source,  circu¬ 
lation,  and  sediment  dynamic  history.  The  trace  metal  iron  is  enriched  in 
riverborne  particles  relative  to  marine  particles ,  and  occurs  at  5-7  percent 
by  dry  weight  (Eaton,  Grant,  and  Gross  1980).  Dissolved  iron  becomes  bound  to 
particulates  at  low  salinities,  and  total  iron  remains  conservative  through 
estuarine  systems.  Ratios  of  manganese,  zinc,  and  other  trace  metals  have 

Table  3 

Trace  Metals  on  TSM 


Tidal 

Concentration,  mg/g 

drv  weight 

Range 

Phase 

Fe 

Mn 

Zn 

Cu 

Ni 

1 

Flood  (F) 

61.4 

1.06 

0.142 

0.092 

0.138 

Ebb  (E) 

60.0 

1.02 

0.185 

0.087 

0.133 

(F  -  E) 

1.4 

0.04 

-0.043 

0.005 

0.005 

2 

Flood 

41.0 

0.816 

— 

0.083 

0.101 

Ebb 

44.8 

0.872 

0.178 

0.076 

0.106 

(F  -  E) 

-3.8 

-0.056 

- 

0.007 

0.005 

3 

Flood 

44.4 

0.843 

0.152 

0.069 

0.107 

Ebb 

42.5 

0.854 

0.149 

0.067 

0.121 

(F  -  E) 

1.9 

-0.011 

0.003 

0.002 

-0.014 

4 

Flood 

39.8 

0.686 

0.139 

0.057 

0.092 

Ebb 

46.4 

0.828 

0.175 

0.069 

0.107 

(F  -  E) 

-6.6 

-0.142 

-0.036 

-0.008 

-0.015 
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Table  4 

Normalized  Trace  Metals  on  TSM 


Tidal 

Fe 

Trace  Metal 

Ratios 

Range 

mg/g  drv  weight 

Mn/Fe 

Zn/Fe 

Cu/Fe 

Ni/Fe 

1 

Flood 

61.4 

0.0173 

0.0023 

0.0015 

0.0022 

Ebb 

60.0 

0.0170 

0.0031 

0.0015 

0.0022 

2 

Flood 

41.0 

0.0199 

_ 

0.0020 

0.0025 

Ebb 

44.8 

0.0195 

0.0040 

0.0017 

0.0024 

3 

Flood 

44.4 

0.0190 

0.0034 

0.0016 

0.0024 

Ebb 

42 . 5 

C.C201 

0.0035 

0.0016 

0.0028 

4 

Flood 

39.8 

0.0172 

0.0035 

0.0014 

0.0023 

Ebb 

46.4 

0.0178 

0.0038 

0.0015 

0.0023 

been  used  to  examine  the  origins 

of  estuarine 

particulates.  Helz  et  al. 

(1985)  gave  the  following  ratios 

for  zinc  and 

manganese : 

Source 

Mn/Fe 

Zn/Fe 

Continental  crust 

0.017 

0.0012 

Sedimentary  rock 

0.020 

0.0028 

Soils 

0.022 

0.0013 

46.  Dissolved  manganese  is  high  (3-12  mg/f)  in  bed  pore  water,  and 
increases  in  the  flow  where  resuspension  occurs  (Morris,  Bale,  and  Howland 
1982).  Thus  manganese  is  nonconservative  due  to  remobilization  of  dissolved 
manganese  from  bed  sediments.  Manganese  returns  to  bed  sediments  by  the 
deposition  of  suspended  sediments.  The  fraction  of  manganese  in  fine  silts 
and  clay  (<16  n)  has  been  found  to  be  a  constant  for  a  particular  sediment 
source  and  may  vary  between  sources  (deGroot  1964). 

47 .  Only  a  few  inferences  can  be  drawn  from  the  trace  metals  data 
presented  in  Tables  3  and  4.  The  strongest  riverine  influence  was  detected  at 
Range  1,  as  indicated  by  the  highest  iron  concentrations  (6  percent  by  dry 
weight) .  Particulate  iron  was  not  conservative  with  resrect  to  the  fraction 
of  fresh  water  as  suggested  for  other  estuarine  systems.  The  influence  of  an 
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alternate  sediment  source  (non"iverine)  is  indicated  for  Rang  •«  2-4  where 
particulate  iron  was  appreciably  reduced.  The  origin  of  the  alternate  sedi¬ 
ment  source  and  the  extent  of  its  influence  could  not  be  determined  from  the 
data.  Estuarine  and  oceanic  sediments  are  both  possible  sources. 

Near-Bed  Transport  Inferred  from  Bed  Forms 

48.  The  net  direction  of  bottom  sediment  transport  was  determined  in  a 
previous  study  (Rubin  and  McCulloch  1979)  using  observations  of  the  movements 
and  shapes  of  sand  bed  forms.  A  100-kHz  side-scan  sonar  was  towed  along  known 
tracks  to  detect  bed  forms.  Some  track  lines  were  repeated  over  several  tidal 
cycles,  and  some  fixed  points  were  monitored  over  longer  periods.  However, 
the  exact  dates  of  various  surveys  were  not  reported,  and  it  is  not  known 
whether  these  conditions  vary  between  summer  and  winter  conditions,  or  with 
tidal  conditions. 

49.  Transport  directions  wire  inferred  from  the  orientations  of  crests 
and  asymmetry  of  sand  waves.  Sand  ripples  were  too  small  to  detect.  Sand 
waves  covered  about  half  of  the  area  in  depths  greater  than  about  10  m  in 
Central  Bay.  Flat  bed  was  the  second  most  abundant  form.  Flat  beds  were  bare 
rock  areas  near  the  Golden  Gate  and  low-velocity  areas  along  the  eastern  edge 
of  Central  Bay.  Transport  directions  based  on  sand  wave  migration  direction 
were  determined  for  the  areas  of  sand  waves.  Figure  16  shows  a  plot  of  the 
results.  Arrows  indicate  direction  of  dominant  transport  (not  magnitude),  and 
lines  indicate  boundaries  between  ebb-  and  flood-dominated  areas.  The  pat¬ 
terns  indicated  in  Figure  16  are  consistent  with  a  flood-dominated  channel 
area  east  of  the  Golden  Gate  flanked  by  ebb-dominated  areas  to  the  north  and 
south.  The  Golden  Gate  appears  to  contain  a  convergence  zone  in  the  south  and 
a  divergence  zone  in  the  north. 

50.  The  Alcatraz  disposal  area  is  indicated  as  ebb  dominated.  Previous 
corings  at  the  Alcatraz  disposal  site  indicated  a  westward  encroachment  of 
sand  in  the  sedimentation  record.*  As  presented  later,  various  current 
velocity  studies  confirm  the  Alcatraz  disposal  area  to  have  ebb-dominant 
maximum  currents. 


*  R.  Farmer.  1985.  "Geotechnical  Investigation,  Alcatraz  Disposal  Site,  San 
Francisco  Bay,  California, 11  Draft  Report,  US  Army  Engineer  District,  San 
Francisco,  San  Francisco,  CA. 
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Figure  16.  Direction  of  near-bottom  current  maximums 
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Near-Bed  Current  Dominance 


51.  Uinzler  and  Kelly  (WK)  (1985),  Science  Applications  International 
Corporation  (SA)  (1987a  and  1987b),  and  the  US  Geological  Survey  (USGS)  (Cheng 
and  Gartner  1984)  have  previously  reported  on  near-bed  currents  based  on  field 
programs.  WK  and  SA  studies  were  commissioned  by  the  San  Francisco  District, 
and  their  measurements  were  carried  out  near  the  Alcatraz  disposal  site.  USGS 
made  measurements  throughout  the  San  Francisco  Bay  system,  including  Central 
Bay. 

Alcatraz  disposal  site  area 

52.  WK  performed  current  measurements  betwem  5  and  11  July  1985.  Four 
current  meter  arrays,  designated  WK-A,  WK-C,  WK-D,  and  WK-E  were  deployed 
within  a  1,500-ft  radius  of  the  Alcatraz  disposal  site.  Sta  WK-A  was  located 
at  the  center  of  the  site,  and  sta  WK-C  through  WK-E  were  located  at  the  east, 
south,  and  west  compass  points  of  the  site  perimeter,  respectively.  A  fifth 
array  (WK-F)  was  located  900  ft  westward  from  the  west  edge  of  the  disposal 
site.  At  each  array  Aanderaa  RC-4  current  meters  were  positioned  at  a  4-m 
depth,  middepth,  and  1.5m  off  the  bottom.  At  sta  WK-A  an  additional  instru¬ 
ment  system  consisting  of  a  wave  recorder,  electromagnetic  current  meter,  and 
transmissometer  was  deployed  0.75  m  above  the  bottom. 

53.  Table  5  presents  means  and  5  percent  exceedences  for  flood  and  ebb 
current  speeds  measured  by  the  WK  current  meters.  In  all  cases  the  ebb  cur¬ 
rent  speeds  are  greater  than  flood  speeds. 

54.  SA  performed  current  measurements  near  the  Alcatraz  disposal  site 
between  4  July  and  10  December  1986.  General  Oceanics  6011  MKI  current  meters 
(Niskin  standard  wing-type  recording  tilt)  were  deployed  1.5  m  above  the  bed 
at  three  stations:  SA-E  (east);  SA-W  (west);  and  SA-S  (south).  The  same  type 
meters  were  deployed  at  2-m  depth,  12-m  depth,  and  1.5-m  above  the  bed  at 

sta  SA-3M.  A  Sea  Data  electromagnetic  current  meter  and  Sea  Tech  transmissom¬ 
eter  were  deployed  at  sta  SA-tripod  1.5  m  above  the  bed.  See  Figure  17  for 
current  meter  locations. 

55.  Table  6  presents  5  percent  exceedence  maximum  currents.  Ebb  cur¬ 
rent  speeds  are  generally  greater  than  corresponding  flood  current  speeds, 
although  at  SA-S  and  SA-3M  (bottom) ,  maximum  ebb  and  flood  current  speeds  are 
similar.  Sta  SA-W  was  located  west  of  the  disposal  site  and  recorded 
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Table  5 

Currents  Near  Alcatraz  Disposal  Site 
14-21  July  1985  (Winzler  and  Kelly  1985) 


Station 


(Station  Depth) 
m 

Position 

Meter 

Depth 

m 

WK-A 

Top 

4 

(13.7) 

Middle 

6 

Bottom 

12 

Tripod 

13 

WK-C 

Top 

4 

(12.5) 

Middle 

6 

Bottom 

11 

WK-D 

Top 

4 

(21.0) 

Middle 

10 

Bottom 

20 

WK-E 

Top 

4 

(16.5) 

Middle 

8 

Bottom 

15 

WK-F 

Bottom 

23 

(24.7) 

Current  Sneed,  cm/sec 


Flood 

Ebb 

5%  Exceeded 

Mean 

5%  Exceeded 

Mean 

80 

51 

114 

62 

77 

50 

104 

59 

67 

43 

86 

48 

62 

40 

83 

44 

81 

51 

123 

64 

95 

56 

130 

71 

70 

43 

90 

50 

90 

57 

125 

65 

93 

63 

129 

74 

54 

37 

71 

40 

64 

33 

73 

36 

87 

57 

116 

67 

60 

38 

68 

40 

36 

20 

40 

21 
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Table  6 

Current  Near  Alcatraz  Disposal  Site 
4  Julv-10  December  1986  (SA  1987a) 


Current  Speed,  cm/sec,  at  5%  Exceedence 


Flood _  _ Ebb 


Sta 

Position 

7-4 

7-28 

7-31 

9-4 

9-9 

10-30 

11- 7 

12- 10 

7-4 

7-28 

7-31 

9-4 

9-9 

10-30 

11- 

12- 

SA-E 

Bottom 

59 

* 

60 

59 

70 

* 

64 

69 

SA-W 

Bottom 

97 

105 

88 

95 

80 

90 

74 

75 

SA-S 

Bottom 

84 

77 

84 

86 

93 

86 

86 

85 

SA-3M 

Surface 

107 

* 

91 

* 

129 

* 

115 

* 

SA-3M 

Middepth 

89 

* 

73 

* 

102 

* 

84 

* 

SA-3M 

Bottom 

63 

* 

59 

* 

64 

* 

56 

* 

SA-tripod 

Bottom 

62 

66 

44 

49 

75 

69 

48 

54 

Note:  Top  figure  in  each  column  is  the  deployment  date;  bottom  figure  is  the 
retrieval  date. 

*  No  data  available. 
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greater  flood  currents  than  ebb,  consistent  with  the  pattern  suggested  in 
Figure  16 . 

Central  Bay  tidal  dominance 

56.  Tides  in  Central  Bay  are  a  nearly  equal  mixture  of  diurnal  and 
semidiurnal  harmonic  tidal  constituents.  Tides  have  two  high  waters  and  two 
low  waters  like  semidiurnal  tides,  but  with  unequal  amplitudes.  Higher  high 
water  (hhw)  normally  precedes  lower  low  water  (llw).  Tidal  currents  following 
this  pattern  would  be  expected  to  be  strongest  during  this  hhw-llw  ebb  tidal 
phase . 

57.  Tidal  current  observations  from  a  large  data  set  collected  in  1979- 
1980  (Cheng  and  Gartner  1984)  were  examined  to  determine  tidal  dominance. 
Evaluating  maximum  current  speeds  by  harmonic  analysis  technique  requires  a 
minimum  record  length  of  14  days  and  preferably  28  days  or  longer.  Daily 
tidal  dominance  can  change  during  a  lunar  month. 

58.  Tidal  dominance  was  estimated  by  computing  the  relative  phase 
between  the  M2  (the  principal  harmonic  constituent)  and  the  M4  shallow-water 
harmonic  constituent  (similar  to  Friedrichs  and  Aubrey  1988)  .  Using  twice  the 
M2  phase  less  the  M<,  phase,  symmetric  tidal  currents  were  defined  as  having 
relative  phases  of  about  90  or  270  deg  while  ebb  dominated  were  between  90  and 
270  deg  and  flood  dominated  between  270  deg  clockwise  through  90  deg.  Current 
data  plots  were  spot  checked  to  confirm  these  relationships. 

59.  Bottom  current  results  are  shown  in  Figure  18  for  long  time  series 
(greater  than  28  days).  Results  were  not  always  consistent  with  other  analy¬ 
ses.  For  example,  the  area  south  of  Alcatraz  Island  was  indicated  as  flood 
dominant,  while  other  analyses  in  this  general  area  indicated  ebb  dominance. 

Hydraulic  Residence  Time 

60.  This  section  evaluates  hydraulic  residence  time  characteristics  of 
the  Central  Bay  area.  Hydraulic  residence  time  is  the  average  time  for  a 
parcel  of  water  initially  residing  in  a  certain  area  or  volume  of  concern  to 
leave  permanently.  Previous  physical  model  test  results  were  reanalyzed  and 
used  for  this  purpose.  The  purpose  of  reviewing  hydraulics  residence  time 
characteristics  is  primarily  to  examine  flow  and  mixing  regimes  in  a  descrip¬ 
tive  way,  and  to  gain  insight  into  sediment  dispersal  in  the  Central  Bay  area. 

61.  When  injected  into  model  or  prototype  estuaries,  dyes  and  other 
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SAN  FRANCISCO 


122°-25'W 


Figure  18.  Tidal  dominance  of 


122°-20'W 

bottom  currents 


conservative  tracers  follow  the  flow  and  do  not  behave  exactly  like  sediments. 
Sediments  do  not  mix  irreversibly  with  water  as  dyes  do.  They  readily  settle 
and/or  stratify  vertically  during  flow  slacks,  and  can  concentrate  in  areas 
and  in  ways  different  from  dye.  Still,  dye  tracers  are  representative  of 
nonsettling  sediment  wash  loads.  Dye  tracer  tests  can  be  descriptive  of  the 
same  flow  and  mixing  regimes  that  strongly  influence  sediment  transport. 

62.  Dye  dispersion  tests  were  previously  performed  in  the  physical 
hydraulic  model  of  the  SFBM  operated  by  the  San  Francisco  District.  (The  SFBM 
is  described  in  Part  IV.)  The  tests  were  performed  for  the  purpose  of  tracing 
the  movements  of  pollutants  and  gaging  the  effects  of  proposed  tidal  barriers 
(USAED,  San  Francisco,  1963).  A  net  delta  outflow  of  16,000  cfs  and  19-year 
mean  tide  were  used  for  these  tests. 

63.  Dye  loss  and/or  decay  was  accounted  for  in  the  model  tests. 

Pontacyl  Brilliant  Pink-B  dye,  an  acid  form  of  Rhodamine-B  dye,  was  found  to 
be  most  conservative  and  was  used  in  the  tests.  Special  tests  were  conducted 
in  several  areas  of  the  model  to  determine  dye  loss.  Dye  concentration  data 
for  a  given  tidal  cycle,  without  respect  to  location  in  the  hydraulic  model  or 
initial  concentration  of  the  dye  release,  were  averaged  to  construct  correc¬ 
tion  curves.  Individual  dye  loss  curves  varied  by  less  than  5  percent  during 
30  tidal  cycles.  After  40  cycles,  about  83  percent  of  the  initial  dye  was 
retained  in  the  model  according  to  the  average  correction  curve. 

64.  For  each  test,  a  solution  of  dye  was  mixed  with  seawater  to  the 
tidal  average  density  of  flow  at  the  injection  point.  Dye  pulses  were  in¬ 
jected  over  a  lunar  day  at  uniform  rates  beginning  about  12  hr  before  local 
hhw.  Injections  were  made  at  the  surface  but  were  seen  to  mix  vertically 
throughout  the  depth.  Samples  were  taken  at  times  of  local  hhw  slack  current. 
Between  1,500  and  3,000  samples  were  generally  collected  during  each  test. 
Three  data  sets  were  used  in  analyses  presented  here  with  injections  as 
follows : 


Ini ection 
I 
R 
H 


_ Shoreside  Location _ 

Passenger  Pier,  San  Francisco 
Point  Richmond,  Richmond 
Islais  Creek,  South  Bay 


65.  Dye  test  results  were  used  to  calculate  residence  times  n  and 
transition  probabilities  Pi:j  .  A  compartmental  analysis  was  used. 
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Central  Bay  and  adjoining  areas  were  segmented  into  five  areas  as  shown  in 
Figure  19 ,  and  defined  as  follows : 

Segment  Description  Area,  so  ft 


1 

Lower  San  Pablo  Bay 

0.91E9 

2 

Upper  Central  Bay 

1.76E9 

3 

Offshore  bar  area 

1.89E9 

4 

Lower  Central  Bay 

1.66E9 

5 

Upper  South  Bay 

1.57E9 

66.  Three  phases  of  dispersion  can  be  identified  in  the  physical  model 
dye  tests.  During  the  short  initial  phase  of  one  to  two  tidal  cycles,  the  dye 
was  advected  in  parcels  at  high  concentration  by  tidal  currents,  usually  into 
all  five  segments.  During  the  second  phase,  local  longitudinal  concentration 
gradients  were  large,  and  dye  concentrations  fell  rapidly  especially  near  the 
injection  segment.  In  the  third  phase,  after  5  to  10  tidal  cycles,  concentra¬ 
tions  were  more  uniform  and  dye  concentrations  fell  slowly. 

67.  In  a  one-dimensional  flow-through  system,  the  local  flushing  rate 
Uj  can  be  defined  by  tracer  experiment  as: 


U.  = 


(2) 


where  M  is  the  mass  of  tracer  injected  upstream,  t  is  time,  and  Cj  is 
the  average  concentration  (mass  per  unit  volume)  in  segment  j .  Then: 


7T  = 


(3) 


where  V^-  is  the  volume  of  the  segment.  In  this  method,  the  integral  dye 
concentration  history  at  a  point  is  used  to  define  a  flushing  rate  for  a  seg¬ 
ment,  implying  a  compartment  analysis.  An  alternate  approach  might  be  to 
calculate  dispersion  coefficients  from  the  second  moments  of  the  dye  concen¬ 
tration  histories.  However,  second  moments  are  sensitive  to  the  distribution 
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tails,  and  the  test  data  are  of  insufficient  length  and  temporal  resolution  to 
perform  such  computations  accurately. 

68.  Some  additional  considerations  are  required  for  the  analysis  of 
San  Francisco  Bay  with  its  unsteady  flow  and  its  multiple  connections  or 
branches.  The  local  flushing  rate  under  these  conditions  (Zvirin  and  Shinnar 
1976)  can  be  expressed  as: 


It  = 


w 


where  P^  is  the  probability  that  a  particle  launched  in  segment  i  would 
pass  through  segment  j  at  least  once,  and  Mi  is  the  mass  launched  at  i  . 
Pi:j  can  be  evaluated  from  the  tracer  tests.  In  each  segment  j  ,  transition 
probabilities  were  calculated  as: 


Pij(t)  = 


Cj(t)  V3(t) 
Mi 


(5) 


where  pi:j(t)  is  the  fraction  of  particles  residing  in  segment  j  and  re¬ 
leased  from  segment  i  at  time  t  .  Then: 


P 


id 


J"  pi5(t)  dt 

JT  PM<C)  dt 


(6) 


Theoretically,  Px:j  for  any  segment  between  the  injection  segment  and  a  sin¬ 
gle  outlet  would  be  unity,  and  P^  =  1;  however,  the  interpretation  of  7r 
for  a  branched,  unsteady  system  is  more  complicated  than  that  for  a  linear 
flow-through  system,  as  will  be  discussed  later. 

69.  Concentration  and  mass  of  dye  for  the  five  segments  were  estimated 
for  each  of  the  model  data  collection  periods  1,  3,  5,  10,  20,  30,  and  40 
tidal  cycles  after  injections  I  and  R,  and  1,  3,  5,  10,  and  20  cycles  after 
injection  H  (a  shorter  test).  Mass  values  for  tidal  cycles  1  and  3  were  used 
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to  extrapolate  back  to  time  zero  and  estimate  the  initial  mass  injected. 

70.  Dye  test  results  from  the  SFBM  generally  indicated  that  the  main 
portion  of  Central  Bay  was  a  well-flushed  section  of  San  Francisco  Bay,  and 
that  dye  releases  resulted  in  low  concentrations  with  respect  to  the  amount  of 


released. 

Results  of  the  n 

and  Pjj 

analyses 

were  as 

follows : 

7r,  Tidal  Cycles  for 
Iniection  Test 

Pij  for 
Iniection 

Test 

Segment 

I  H 

R 

I 

H 

R 

1 

0.0  0.1 

1.4 

0.1 

0.3 

0.6 

2 

0.7  0.7 

4.1* 

0.6 

0.7 

1.0* 

3 

2.5  1.4 

1.7 

1.2 

1.0 

0.6 

4 

1.8*  1.9 

6.5 

1.0* 

1.2 

1.3 

5 

0.6  1.4* 

0.6 

0.5 

1.0* 

0.4 

*  Indicates  injection  segment. 

71.  In  a  simple  system  (not  San  Francisco  Bay)  the  value  of  ir  should 
be  independent  of  injection  location,  and  hence  should  be  consistent  between 
tests.  Some  of  the  variation  in  tt  results  can  be  explained  by  local  effects 
near  injection  points  R  and  H.  Following  injection  R,  high  concentrations 
persisted  in  the  shallow  areas  near  Richmond,  and  were  subsequently  trans¬ 
ported  into  the  shallow  area  off  Berkeley.  Dye  was  slow  to  disperse  into  the 
main  portion  of  the  flow  from  these  shallow  locations.  The  initial  dye  mass, 
estimated  as  described  earlier,  was  only  about  half  that  of  other  tests,  indi¬ 
cating  that  some  of  the  dye  may  have  stayed  near  the  injection  point  and  was 
not  immediately  sampled.  Residence  times  for  segments  2  and  4  were  much 
larger  for  this  test,  and  reflected  the  slower  flushing  rates  of  the  shallow 
area  of  these  segments.  The  injection  at  H  similarly  produced  highest  concen¬ 
trations  in  shallow-water  areas,  in  this  case  near  Sierra  Point  in  South  Bay, 
and  the  greatest  7r  for  the  injection  segment  (5). 

72.  These  shallow-water  effects  have  implications  with  respect  to  the 
dispersion  of  wave-resuspended  sediments  from  those  areas.  Flushing  of  dye 
from  the  shallow  area  off  Berkeley  in  the  SFBM  was  particularly  slow.  There¬ 
fore,  resuspended  sediments  would  have  the  opportunity  to  redeposit  before 
being  mixed  with  the  main  portions  of  the  flow  and  transported  away  from  this 
area.  This  applies  to  shallows  on  the  west  side  of  upper  South  Bay  to  a 
lesser  extent.  The  slow  flushing  of  the  shallow-water  areas,  if 
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representative  of  the  prototype,  would  reduce  the  influence  of  these  areas  and 
the  processes  affecting  them  (i.e.,  wind  waves)  on  overall  transport  of 
suspended  sediments. 

73.  The  extremely  short  values  of  tt  for  lower  San  Pablo  Bay  (seg¬ 
ment  1)  are  unrealistic,  a  limitation  of  the  analysis  pointed  out  by  Zvirin 
and  Shinner  (1976).  Segment  1  was  poorly  connected  to  the  other  regions  ac¬ 
cording  to  values  of  P^  .  Only  a  small  amount  of  mass  reached  segment  1, 
which  had  strong  local  mixing. 

74.  As  discussed  earlier,  Pi:j  is  dependent  on  injection  location,  has 
a  theoretical  maximum  of  1.0,  and  should  be  1.0  in  the  injection  segment  and 
in  all  segments  downstream  from  the  injection  segment.  The  calculated  Pi;j 
indicates  a  strong  connection  between  segments  3-5  for  injection  H.  Pi;j  has 
a  theoretical  maximum  of  1.0,  and  the  occurrence  of  values  greater  than  1.0 
indicate  that  model  test  results  were  probably  accurate  only  to  within 

±20  percent  or  so. 

75.  The  injection  at  I  produced  the  most  immediate  dispersion.  Maximum 
concentrations  were  in  the  main  portions  of  the  flow.  Results  of  v  for 
injection  H  were  similar  to  I  except  for  segment  5,  which,  as  explained 
earlier,  was  affected  by  low  shallow-water  flushing. 

76.  Local  residence  times  are  descriptive  of  flushing  conditions,  and 
can  also  be  used  to  calculate  equilibrium  concentrations  resulting  from  a 
continuous  discharge.  The  results  for  the  injection  at  I  were  used  to 

A 

estimate  the  buildup  of  turbidity  C  caused  by  a  constant  release  of  mate¬ 
rial.  The  previously  mentioned  qualifications  apply.  A  constant  release  rate 
R0  would  roughly  correspond  to  the  continuous  release  of  a  conservative, 
passive  tracer  at  the  Alcatraz  disposal  site.  Then  the  release  rate  into  any 
segment  is 


Rj  ■  Ko 


(7) 


where  a ^  is  a  distribution  parameter,  and 
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Chen  at  steady  state, 


Cs  =  JTj 

vj 


(8) 


Results  for  a  10,000-metric-ton  release  per  tidal  cycle  were  as  follows: 


Segment 

_ j _ 

1 

2 

3 

4 

5 


C 

0.0 

0.7 

2.3 

2.0 

0.9 


Caution  should  be  exercised  when  interpreting  these  results.  As  described 
previously,  dye  is  a  poor  analog  for  sediment.  The  physical  model  from  which 
these  data  were  developed  was  never  verified  to  residence  times,  and  scale 
distortion  effects  may  have  affected  dye  test  results.  The  only  cross-check 
on  the  estimated  residence  times  on  hand  is  the  prototype  hydraulic  residence 
time  estimated  by  dividing  segment  volumes  by  tidal  prism.  Tidal  prisms  are 
available  only  at  the  three  boundaries  of  Central  Bay,  and  segments  2  and  4 
had  to  be  combined  for  this  analysis.  Using  a  mean  tide,  the  hydraulic 
residence  time  of  Central  Bay  (segments  2  and  4)  was  1.6  tidal  cycles. 
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PART  IV:  DESCRIPTION  OF  THE  MODELS 


Model  Overview  and  Procedures 


77.  The  prediction  of  TSM  concentration  in  Can  f rancisco  Bay  was  accom¬ 
plished  with  mathematical  models.  A  hydrodyna..J<-  mode!  generated  the  time- 
varying  currents  and  water  elevations  at  compute t lonal  nodes  in  a  numerical 
mesh  representing  San  Francisco  Bay.  These  hydrodynamics  were  used  in  a  sedi¬ 
ment  transport  model  to  solve  the  convection-diffuse....!  equation  with  appropri¬ 
ate  sink/source  terms  for  TSM.  The  two  numerical  models,  RMA-2V  (A  Two- 
Dimensional  Model  for  Free  Surface  Flows)  and  STUDH  (Sediment  Transport  in 
Unsteady  2-Dimensional  Flows,  Horizontal  Plane),  are  included  in  the  TABS 
modeling  system,  which  is  supported  by  the  US  Army  Corps  of  Engineers  (Thomas 
and  McAnally  1985).  Data  from  the  SFBM  were  the  principal  source  of  verifica¬ 
tion  currents  and  water  levels  for  the  hydrodynamic  model. 

The  San  Francisco  Bay-Delta  Physical-Hydraulic  Model 

78.  The  SFBM  is  a  fixed-boundary,  distorted-scale,  tidal  hydraulic 
model.  This  physical  model  provides  a  means  of  reproducing,  on  a  manageable 
scale,  the  three-dimensional  hydrodynamic  phenomena  that  occur  throughout  the 
large  and  complex  estuarine  system.  Similarity  between  model  and  prototype  is 
established  using  laws  of  similitude  expressed  by  a  set  of  dimensionless  num¬ 
bers.  These  may  include  the  Froude,  Reynolds,  Weber,  and  Cauchy  numbers, 
which  are  the  ratios  of  inertial  forces  to  the  forces  of  gravity,  viscosity, 
surface  tension,  and  elasticity,  respectively.  Perfect  similitude  would  re¬ 
quire  that  each  of  these  numbers  be  the  same  for  model  and  prototype.  Physi¬ 
cal  hydraulic  models  cannot  be  designed  to  provide  similitude  with  respect  to 
all  forces,  but  are  designed  to  simulate  the  dominant  forces  affecting  the 
conditions  being  modeled.  This  requires  different  horizontal  and  vertical 
scales,  and  therefore  a  distorted-scale  model.  The  scales  were  selected  so 
that  the  remaining  forces  are  negligible  or  do  not  cause  fundamental  dissimi¬ 
larities  in  the  processes  of  interest. 

79.  In  the  San  Francisco  Bay  complex,  the  depth,  surface  slope,  and 
other  features  of  tidal  flow  are  controlled  by  the  joint  effects  of  inertial 
and  gravitational  forces.  The  major  hydraulic  quantities  vary  according  to 
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the  Froude  number.  For  the  SFBM,  different  model  scales  were  chosen  for 
horizontal  and  vertical  dimensions  to  reduce  construction  and  operation  time 
and  cost  as  well  as  reproduce  flows  faithfully  and  permit  measurements  of  the 
desired  parameters  with  satisfactory  accuracy.  The  effect  of  distortion  on 
boundary  roughness  was  reproduced  during  verification  tests  by  the  use  of 
metal  strips  embedded  in  the  model.  These  strips  were  adjusted  to  obtain 
reproduction  of  the  bay's  tides,  currents,  and  salinities.  The  selected 
scaling  ratios  of  model  to  prototype  are  as  follows: 

Scale  Relations 
Characteristic  Model : Prototype 

Horizontal  length  1:1,000 

Vertical  length  1:100 

Time  1:100 

Velocity  (horizontal)  1:10 

Discharge  1:1,000,000 

Salinity  1:1 

80.  The  SFBM  occupies  an  area  of  about  1  acre  and  is  completely  en¬ 
closed  in  a  128,500-sq-ft  shelter  to  protect  it  from  the  weather  and  thus 
permit  uninterrupted  operation.  The  limits  of  the  model,  shown  in  Figure  20, 
encompass  a  portion  of  the  Pacific  Ocean  extending  17  miles  west  of  the  Golden 
Gate,  San  Francisco  Bay,  San  Pablo  Bay,  Suisun  Bay,  and  all  of  the  Sacramento- 
San  Joaquin  Delta  east  of  Suisun  Bay  to  the  cities  of  Sacramento  on  the  north, 
Stockton  on  the  east,  and  Tracy  to  the  south. 

81.  The  model  is  equipped  with  the  necessary  appurtenances  including 
freshwater  and  saltwater  supply  sumps ,  primary  saltwater  inflow  and  mixing 
pumps,  computer-controlled  tide  generator  capable  of  reproducing  repetitive 
ocean  tides  by  regulating  two  motorized  gate  valves,  precise  flowmeters  to 
regulate  major  river  inflows  to  the  delta  plus  other  inflow  and  withdrawal 
facilities,  major  pumping  plant  exports,  and  electrically  operated  one-way 
head  gates.  The  model  represents  the  state  of  the  art  in  model  operation, 
real-time  data  acquisition  using  miniature  sensors,  data  analysis,  and 
information  storage/retrieval  capability. 


The  TABS  Numerical  Modeling  System 


82.  The  TABS  is  a  modular  system  composed  of  distinct  computer  programs 
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Figure  20.  Limits  of  SFBM  (from  USAED,  San  Francisco,  1981) 


linked  together  by  preprocessors  and  postprocessors.  Each  of  the  major  com¬ 
puter  programs  solves  a  particular  type  of  problem:  hydrodynamics  (RMA-2V) , 
sediment  transport  (STUDI1) ,  or  constituent  transport  (RMA-4) .  These  programs 
employ  the  finite  element  method  to  solve  the  governing  equations.  A  brief 
description  of  RMA-2V  and  STUD1I  appears  in  Appendix  B. 

Numerical-Hvdrodvnamic  Model 

83.  RMA-2V  is  a  finite  element  solution  of  the  Reynolds  form  of  the 
Navier-Stokes  equations  for  turbulent  flows.  Friction  is  calculated  with 
Manning's  equation,  and  eddy  viscosity  coefficients  are  used  to  define  turbu¬ 
lent  exchange  characteristics.  A  velocity  form  of  the  basic  equation  is  used 
with  side  boundaries  treated  as  either  slip  (parallel  flow)  or  static  (zero 
flow) .  The  model  recognizes  computationally  wet  or  dry  elements  and  corrects 
the  mesh  accordingly.  Boundary  conditions  may  be  water-surface  elevations, 
velocities,  or  discharges  and  may  occur  inside  the  mesh  as  well  as  along  the 
outer  boundaries. 


Numerical  Mesh 


Model  limits 

84.  A  comprehensive  mesh  of  the  San  Francisco  Bay  system  was  developed 
using  the  original  global  mesh  from  the  previous  numerical  study  of  the 
Alcatraz  disposal  site  reported  in  Pankow  (1988).  The  new  mesh  includes  ex¬ 
pansions  to  the  Pankow  (1988)  global  mesh  through  the  two-dimensional  repre¬ 
sentation  of  Suisun  Bay,  a  larger  ocean  boundary  in  the  Gulf  of  the  Faralones 
(the  Pacific  Ocean),  and  the  one-dimensional  representation  of  the  Napa  River, 
Suisun  Slough,  Montezuma  Slough,  and  the  Sacramento-San  Joaquin  Delta  (Fig¬ 
ure  21) .  The  numerical  representation  of  the  Pacific  Ocean  encompasses  the 
region  between  Drakes  Bay  to  the  north,  the  Farallon  Islands  to  the  west,  and 
Half  Moon  Bay  to  the  south,  which  are  all  approximately  22  nautical  miles  from 
the  Golden  Gate  Bridge.  The  mesh  consists  of  2,330  elements  with  7,049  nodes. 
Quadrilaterals  and  triangles  compose  the  2,019  two-dimensional  elements  repre¬ 
senting  the  open  bay  and  ocean  waters,  while  311  one-dimensional  elements 
define  the  riverine  portions  of  the  system.  This  mesh  is  referred  to  as 
Mesh  A. 
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ierical  model  representation  of  San  Francisco  Bay  (Mesh  A  and  Mesh  B) 


85.  Because  the  sediment  transport  model  is  presently  limited  to  using 
two-dimensional  elements,  a  second  mesh  (Mesh  B)  was  developed  that  was  iden¬ 
tical  to  the  expanded  comprehensive  mesh  except  for  the  exclusion  of  the 

311  one-dimensional  elements  (Figure  21) .  The  method  of  using  both  meshes  is 
described  in  the  boundary  conditions  section.  A  second  reason  for  using 
Mesh  B  is  the  lack  of  resolution  in  the  one-dimensional  portion  of  Mesh  A  to 
address  adequately  sedimentation  and  TSM  values  in  the  Sacramento-San  Joaquin 
Delta. 

Bathymetry 

86 .  The  bathymetric  data  required  by  the  numerical  models  were  obtained 
from  the  following  NOS  nautical  charts: 


Chart  No. 

Location 

Scale 

Date 

18645 

Gulf  of  the  Faralones 

1:100,000 

1976 

18651 

San  Francisco  Bay-southern  part 

1:40,000 

1978 

18649 

Entrance  to  San  Francisco  Bay 

1:40,000 

1986 

18654 

San  Pablo  Bay 

1:40,000 

1988 

18656 

Suisun  Bay 

1:40,000 

1983 

A  state  coordinate 

(x-coordinate  and  y-coordinate)  and  a  z-value  (bed  eleva- 

tion)  were  assigned 

to  each  node  in  both  Meshes  A  and  B  through 

use  of  a 

digitizer.  The  vertical  datum  was  set  to  mean  lower 

low  water 

(mllw) . 

Time-step 

87.  The  time-step  used  in  the  hydrodynamic  numerical  simulations  is 
identical  to  that  used  in  the  SFBM  for  data  collection.  This  allows  the 
direct  comparison  of  physical  model  and  numerical  model  results.  A  time-step 
of  37.26  min  (0.621  hr)  was  used,  which  provides  40  equal  steps  in  the 
24.84-hr  lunar  day. 

Boundary  conditions 

88.  If  no  boundary  condition  is  specified  at  a  node,  RMA-2V  computes  an 
x-  and  y-component  of  velocity  and  a  water  depth.  For  boundaries  along 
land/water  interfaces,  slip  flow  parallel  to  the  boundary  is  specified.  Along 
the  Pacific  Ocean  outer  boundary,  a  time-dependent  head  specification  (water- 
surface  elevation)  was  assigned  to  all  boundary  nodes  to  represent  tidal  ele¬ 
vation  fluctuations.  In  the  Sacramento-San  Joaquin  Delta,  velocity  specifica¬ 
tions  were  made  at  the  Mesh  A  headwaters  of  the  San  Joaquin  and  Sacramento 
Rivers  and  withdrawals  of  water  were  specified  at  pumping  locations  of  the 
Contra  Costa  Canal  and  the  US  Bureau  of  Reclamation  Tracy  Pumping  Plant.  All 
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other  riverine  freshwater  inflows  were  considered  sufficiently  small  to  be 
assumed  to  be  zero,  which  is  a  reasonable  assumption  for  the  low  inflow  condi¬ 
tions  of  this  study.  In  Mesh  A,  the  Napa  River,  Suisun  Slough,  Sacramento 
River,  and  San  Joaquin  River  one-dimensional  reaches  were  extended  suffi¬ 
ciently  to  include  a  majority  of  the  tidal  prism  associated  with  each  riverine 
system.  While  none  of  these  one-dimensional  riverine  reaches  was  extended 
upstream  beyond  the  influence  of  the  tide,  the  reaches  were  extended  suffi¬ 
ciently  so  that  the  tide  had  appreciably  damped  and  minor  spurious  tidal  re¬ 
flections  induced  by  the  boundaries  had  been  sufficiently  reduced  or  removed 
from  all  areas  of  interest  to  this  study.  A  final  boundary  condition  in 
Mesh  A  was  the  boundary  in  extreme  South  Bay  at  the  Dumbarton  Bridge.  Data 
obtained  from  the  SFBM  were  used  to  supply  the  velocity  boundary  conditions  at 
this  location. 

89.  Mesh  B  used  the  same  time-dependent  water  level  specification  for 
boundary  nodes  in  the  Pacific  Ocean  and  the  same  time-dependent  velocity 
specification  for  boundary  nodes  in  extreme  South  Bay  as  for  Mesh  A.  However, 
with  the  deletion  of  one-dimensional  elements  from  Mesh  B,  new  boundary  spec¬ 
ification  locations  occurred  at  the  termination  of  the  two-dimensional  ele¬ 
ments  in  Mesh  B  at  the  Napa  River,  Suisun  Slough,  Montezuma  Slough,  and 
eastern  Suisun  Bay  (Figure  21) .  The  velocity  specifications  at  these  loca¬ 
tions  were  determined  through  the  laborious  but  necessary  process  of  extract¬ 
ing  them  from  the  results  generated  from  the  execution  of  RMA-2V  for  Mesh  A. 

Numerical-Sediment  Transport  Model 

90.  The  sediment  transport  model  in  the  TABS  system  is  STUDH,  which 
solves  the  convection-diffusion  transport  equation  with  bed  source  terms. 

These  terms  are  structured  for  either  sand  or  cohesive  sediments.  The  Ackers- 
White  (1973)  procedure  is  used  to  calculate  a  sediment  transport  potential  for 
the  sands  from  which  the  actual  transport  is  calculated  based  on  the  avail¬ 
ability.  Clay  erosion  is  treated  as  described  in  Ariathurai,  MacArthur,  and 
Krone  (1977),  and  the  deposition  of  clay  uses  Krone's  (1962)  equations. 
Deposited  material  forms  layers,  and  bookkeeping  within  the  STUDH  code  allows 
up  to  10  layers  at  each  node  for  maintaining  separate  material  types, 
deposited  thickness,  density,  and  age.  The  identical  Mesh  B  used  by  RMA-2V 
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was  used  by  STUDH.  Boundary  conditions  were  specified  as  a  TSM  concentration, 
which  may  be  time-varying. 

Time-s tep 

91.  STUDH  was  operated  with  a  time-step  equal  to  one-half  the  value 
used  in  RMA-2V;  that  is,  80  time-steps  of  18.63  min  each  represented  a  lunar 
day.  The  hydrodynamic  conditions  used  in  STUDH  were  obtained  by  linearly 
interpolating  RMA-2V  results  to  generate  the  hydrodynamic  conditions  for  each 
STUDH  time-step. 

Boundary  conditions 

92.  At  all  boundary  locations  specified  in  RMA-2V,  a  STUDH  TSM  boundary 
concentration  was  supplied.  While  in  actuality  this  boundary  is  undoubtedly 
time  and  spatially  varying  with  tidal  fluctuations  and  other  factors,  insuffi¬ 
cient  data  existed  at  boundary  locations  to  allow  elaborate  representation. 
Instead,  a  single  representative  value  appropriate  for  each  boundary  location, 
e.g.,  Pacific  Ocean  or  Napa  River,  was  determined  and  used  throughout  the 
application  of  STUDH. 
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PART  V:  MODEL  VERIFICATION  TO  SUMMER  CONDITION 


Physical  Model  Verification 

93.  The  verification  report  of  the  San  Francisco  physical  model  (USAED, 
San  Francisco,  1963)  indicated  that  the  model  could  reproduce  a  prototype  mean 
tide  (21-22  September  1956)  with  good  results.  Maximum  error  in  elevation  in 
model  reproduction  at  23  gages  throughout  the  bay  was  0.4  ft  prototype,  with 
the  model  amplitudes  being  slightly  less  than  those  of  the  prototype.  At  most 
stations  the  error  was  only  0.1  ft  prototype.  Water  levels  were  found  to 
repeat  cycle  to  cycle  and  test  to  test.  The  standard  deviation  of  water 
levels  was  found  to  be  0.04  to  0.08  ft  prototype.  Forty  percent  or  more  of 
the  standard  deviation  was  attributed  to  the  tide  generator  (USAED, 

San  Francisco,  1984). 

94.  The  verification  procedure  indicated  that  local  model  velocities 
may  vary  up  to  20  percent  from  those  in  the  prototype,  the  model  velocities 
tending  to  be  generally  higher  than  those  in  the  prototype.  Additional  proto¬ 
type  to  model  comparisons  made  in  1976  (USAED,  San  Francisco,  1976)  were 
similar . 


Hydrodynamic  Model  Verification 

95.  Verification  of  the  RMA-2V  model  employed  in  the  previous  dredged 
material  disposal  study,  which  used  a  global  mesh  of  San  Francisco  Bay,  is 
reported  in  Pankow  (1988).  The  expanded  mesh  (Mesh  A)  for  the  present  study, 
including  the  one-dimensional  representation  in  important  riverine  areas,  was 
verified  to  measured  water  levels  and  velocities  from  the  SFBM  for  mean  and 
spring  tide  conditions.  In  addition,  Mesh  A  was  verified  to  water  levels 
generated  from  a  least  squares  harmonic  analysis  of  prototype  tide  gage  loca¬ 
tions  in  the  San  Francisco  Bay  system.  The  tide  gages  were  operated  and  har¬ 
monic  analysis  performed  as  the  result  of  a  comprehensive  tidal  survey  of  the 
bay  system  conducted  jointly  by  the  USGS  and  NOS/National  Oceanic  and  Atmo¬ 
spheric  Administration  (N0AA)  during  1979  and  1980  (Cheng  and  Gartner  1984; 
Welch,  Gartner,  and  Gill  1985). 

Model  coefficients 

96.  The  process  of  numerical  model  verification  resulted  in  one  final 
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set  of  model  coefficients  valid  for  all  verification  cases  for  Meshes  A  and  B. 
The  model  coefficients  subject  to  adjustment  during  model  verificacion  are 
bottom  roughness  as  represented  by  Manning's  n  value  and  eddy  viscosity 
coefficients.  Based  on  water  depths  and  bottom  conditions,  the  Manning's  n 
values  in  the  Gulf  of  the  Faralones  (Pacific  Ocean)  ranged  from  0.018  to  0.027 
with  the  smaller  values  in  the  deeper  waters.  In  a  similar  manner,  Manning's 
n  values  ranged  from  0.015  in  the  deep  waters  of  San  Francisco  Bay  to  0.025 
in  the  shallow  waters.  In  the  one-dimensional  sections,  slightly  higher 
Manning's  n  values  indicative  of  riverine  conditions  with  a  range  of  0.024 
to  0.035  were  employed.  Eddy  viscosity  values  were  assigned  based  upon  the 
element  size  in  the  mesh.  In  San  Francisco  Bay  the  eddy  viscosity  coeffi¬ 
cients  were  all  set  to  200  lb-sec/ft2;  in  the  one-dimensional  areas  of  the 
Sacramento-San  Joaquin  Delta,  Napa  River,  Suisun  Slough,  and  Montezuma  Slough 
to  1,500  lb-sec/ft2;  and  in  the  Gulf  of  the  Faralones  to  a  range  from  400  to 
4,000  lb-sec/ft2.  Extremely  large  element  areas  in  the  gulf  necessitated  the 
large  eddy  viscosity  values  in  this  area.  The  largest  values  (greater  than 
1,000  lb-sec/ft2)  were  used  in  the  northeast  corner  of  the  mesh  representation 
of  the  Gulf  of  the  Faralones.  The  interaction  of  Coriolis  force  with  a  spec¬ 
ification  of  a  uniform  water  level  along  the  entire  gulf  boundary  resulted  in 
spurious  current  patterns  localized  in  the  elements  in  the  extreme  upper 
(northeast)  gulf  boundary.  This  slight  irregularity  was  controlled  through 
specification  of  large  eddy  viscosities  that  artificially  damped  it.  A 
natural  sloping  water  surface  along  the  ocean  boundary  would  have  accounted 
for  the  influence  of  the  Coriolis  force,  avoiding  this  spurious  current. 

Since  no  data  of  sufficient  detail  exist  to  allow  specification  of  this  slope 
and  since  the  irregularities  were  confined  to  the  extreme  northeast  portion  of 
the  Gulf  of  the  Faralones,  far  removed  from  the  areas  of  interest,  the  damping 
from  eddy  viscosity  was  selected  as  the  most  expedient  manner  to  control  the 
phenomenon . 

Mean  tide  verification.  Mesh  A 

97.  The  numerical  model  (Mesh  A)  was  operated  for  nearly  35  hr  with  the 
19-year  mean  tide  to  produce  a  mixed  tide  with  extremes  of  approximately  0.0, 
4.6,  2.3,  and  5.8  ft  mllw  at  tide  sta  T8,  Golden  Gate  location.  The  ocean 
boundary  conditions  were  specified  to  closely  reproduce  the  tide  at  the  Golden 
Gate  measured  in  the  physical  model.  Net  delta  outflow  in  the  numerical  and 
physical  models  was  specified  at  4,400  cfs.  Numerical  and  physical  model 
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results  for  water  levels  and  velocities  are  compared  in  Plates  1-7.  The  sta¬ 
tion  locations  for  comparisons  are  shown  in  Figure  22.  Comparisons  are  pro¬ 
vided  only  for  selected  stations  that  demonstrate  the  general  trend  of  model 
verification  with  emphasis  on  Central  Bay  stations.  The  complete  set  of  com¬ 
parisons  at  numerous  water  level  and  velocity  stations  is  provided  in  Hauck, 
Heltzel,  and  Teeter  (in  preparation),  the  report  to  the  San  Francico  District 
on  the  salinity  intrusion  and  ship  simulation  study  for  the  proposed  John  F. 
Baldwin  III  Ship  Channel  project. 

98.  Since  the  first  5  to  10  hr  of  numerical  model  results  for  this  mesh 
are  influenced  by  the  model  initial  condition  specification,  only  the  last 
complete  tidal  cycle  (24.84  hr)  of  the  total  35-hr  simulation  is  used  for 
comparison.  Because  the  RMA-2V  code  determines  vertically  averaged  values, 
the  physical  model  data  were  averaged  to  estimate  a  comparable  vertically 
averaged  value.  Generally,  near-bottom  and  near-surface  measurements  were 
taken  in  the  physical  model,  though  only  a  middepth  measurement  was  made  at 
some  shallow  stations  and  for  sta  1  at  the  Golden  Gate  only  a  surface  velocity 
was  measured.  The  mean  velocity  at  sta  1  was  determined  as  85  percent  of  the 
surface  velocity. 

99.  The  water-surface  elevations  for  the  mean  tide  conditions  compare 
closely.  Phasing  of  the  numerical  and  physical  model  is  in  general  agreement, 
though  a  trend  of  RMA-2V  results  preceding  physical  model  measurements  by 
approximately  one-half  hour  is  apparent.  Amplitude  and  tidal  plan  (mean  ele¬ 
vation)  agreement  is  also  favorable.  Only  in  the  Suisun  Bay  area  do  apprecia¬ 
ble  discrepancies  occur  between  the  tidal  planes  in  the  physical  and  numerical 
models.  This  discrepancy  is  the  result  of  the  constant  water  density  limita¬ 
tion  in  RMA-2V,  while  the  physical  model  contains  the  density  effects  (as  in 
the  prototype  system)  of  water  of  varying  salinity  from  least  density  (nearly 
freshwater  conditions)  in  Suisun  Bay  to  greatest  density  (seawater)  in  the 
vicinity  of  the  Golden  Gate.  This  may  be  interpreted  that,  hydrostatically,  a 
higher  level  of  low-density  water  is  required  to  equal  the  pressure  of  higher 
density  water.  The  constant  water  density  in  RMA-2V  manifests  itself  in  an 
underprediction  of  the  tidal  plane  in  the  lower  density  portions  of  the 
system,  i.e.,  Suisun  Bay  and  the  Sacramento-San  Joaquin  Delta. 

100.  Comparison  of  RMA-2V  velocities  to  physical  model  averaged  results 
are  good  considering  the  errors  involved  in  the  process  of  measuring  veloci¬ 
ties  in  the  physical  model  and  errors  resulting  from  representation  of  a 
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Figure  22.  San  Francisco  Bay  area  with  physical  model  station  locations 
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vertical  mean  velocity  by  the  simple  arithmetic  mean  of  a  surface  and  bottom 
measurement.  The  phasing  of  the  numerical  model  velocities  compared  to  phys¬ 
ical  model  velocities  is  good.  At  several  stations,  RMA-2V  results  precede 
the  physical  model  results  by  0.5  to  1  hr,  e.g.,  sta  1A,  IB,  and  1-1B,  while 
at  other  stations  the  phasing  is  very  close,  e.g.,  sta  1-3,  1-4A,  and  1-6A. 
Generally,  the  RMA-2V  velocity  magnitudes  underpredict  the  results  from  the 
physical  model.  However,  as  mentioned  under  physical  model  verification, 
physical  model  velocities  tended  to  be  slightly  higher  than  prototype  veloci¬ 
ties,  so  the  RMA-2V  values  are  considered  to  show  a  good  verification  to 
expected  prototype  velocities. 

101.  The  residual  or  tidally  averaged  currents  for  the  mean  tide  as 
determined  with  RMA-2V  are  depicted  in  Figure  23  for  the  central  San  Francisco 
Bay  area.  With  the  inherent  limitation  that  RMA-2V  results  are  vertically 
averaged,  the  residual  currents  indicate  the  net  direction  and  magnitude  of 
currents  in  an  estuarine  system  where  currents  vary  with  time.  The  residual 
current  pattern  in  Figure  23  is  of  interest  not  only  for  the  predicted  pat¬ 
terns,  but  also  because  of  the  similarity  to  Figure  16,  which  displays  near¬ 
bottom  current  direction  ebb  and  flood  dominance  determined  from  physical 
measurements  and  bed  form  observations.  The  patterns  in  Figures  16  and  23  are 
very  similar,  especially  in  the  areas  to  the  south  and  west  of  Alcatraz 
Island,  around  the  Golden  Gate  Bridge,  in  Raccoon  Strait,  and  south  of 
Richardson  Bay.  An  exception  is  the  bay  area  north  of  Alcatraz  Island  where 
the  numerically  predicted  residual  current  pattern  has  some  notable  discrep- 
encies  from  field  measurements.  However,  RMA-2V  predicted  residual  current 
patterns  are  generally  remarkably  similar  to  near-bottom  patterns  determined 
from  the  physical  measurements,  including  the  locations  of  boundaries  between 
flood  and  ebb  dominance. 

Spring  tide  verification.  Mesh  A 

102.  The  numerical  model  was  operated  in  essentially  the  identical  man¬ 
ner  for  spring  tide  verification  as  for  the  mean  tide  verification.  Again, 
net  delta  outflow  was  specified  as  4,400  cfs  for  numerical  and  physical 
models.  The  numerical  model  was  operated  with  time-varying  water  level  bound¬ 
ary  specifications  in  the  Pacific  Ocean  to  reproduce  the  physical  model  mixed 
tide  with  extremes  of  -1.7,  4.6,  2.3,  and  6.8  ft  mllw  at  tide  sta  T8,  Golden 
Gate . 

103.  The  same  trends  and  characteristics  were  observed  in  the  spring 
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Figure  23.  Numerical  model  predicted  residual  current  pattern  for 

central  San  Francisco  Bay 


tide  verification  as  in  the  mean  tide  verification.  Because  the  verification 
is  similar  and  since  the  mean  tide,  not  the  spring  tide,  was  used  in  the  sedi¬ 
mentation  odeling,  only  a  limited,  representative  number  of  velocity  and 
water  level  verification  stations  are  presented  (Plates  8-10) .  The  complete 
presentation  of  spring  tide  verification  is  presented  in  Hauck,  Heltzel,  and 
Teeter  (in  preparation).  For  the  spring  tide,  tidal  phasing,  amplitude,  and 
plane  are  still  in  good  agreement.  The  half-hour  lag  of  physical  model  phas¬ 
ing  to  numerical  model  phasing  still  occurred,  as  did  the  tidal  plane  discrep¬ 
ancy  in  Suisun  Bay  due  to  the  constant  water  density  assumption  in  RMA-2V. 

The  RMA-2V  velocities  were  in  general  slightly  lower  than  the  vertically  aver¬ 
aged  physical  model  results,  though  numerical  results  were  very  close  or 
greater  than  physical  model  velocities  at  some  stations.  Generally,  the 
phasing  of  the  velocities  showed  reasonable  agreement. 

Harmonic  tide  verification.  Mesh  A 

104.  The  USGS  and  NOS/NOAA  analyzed  tide  (water  level)  records  obtained 
during  their  joint  study  by  least  squares  harmonic  analysis.  The  tide  records 
for  the  various  locations  were  obtained  for  different  lengths  of  time  and  at 
different  times  during  the  1979-80  study  period.  In  an  attempt  to  minimize 
discrepancies  arising  in  the  harmonic  analysis,  the  stations  used  in  verifying 
RMA-2V  were  limited  to  those  with  8  months  or  more  of  water  level  records. 

The  five  major  tidal  constituents  (or  harmonic  constituents)  01,  Kl,  N2,  M2, 
and  S2  were  selected  for  use  in  the  model  verification.  Using  these  five 
major  harmonic  constituents  for  each  station,  the  astronomical  tide  was  pre¬ 


dicted  for  the  period  of  6  and  7  September  1988.  The  astronomical  tide  during 
this  period  was  midway  between  the  extremes  of  spring  and  neap  tides,  there¬ 
fore  representing  an  average  tide . 

105.  The  numerical  model  was  operated  with  no  freshwater  inflow  and 
with  a  tide  that  closely  reproduced  the  harmonically  constructed  tide  at 
sta  T8,  Golden  Gate  (Figure  22).  Comparisons  of  RMA-2V  results  after  model 
spin-up  to  the  harmonically  constructed  tide  at  sta  3,  8,  19,  and  20  are  pro¬ 
vided  in  Plate  11.  Results  of  these  comparisons  are  good.  While  the  RMA-2V 


tide  phasing  still  preceded  the  verification  data  at  sta  3  and  19,  the  differ¬ 


ences  are  not  as  great  as  observed  with  physical  model  results  and  are  gener¬ 


ally  much  less  than  0.5  hr.  At  sta  19,  the  RMA-2V  results  actually  lag  the 


harmonic  tide  by  0.5  hr.  Tidal  amplitudes  from  RMA-2  are  somewhat  too  large, 


and  this  is  most  prominent  at  the  first  high  water.  No  tidal  plane 
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comparisons  could  be  made  since  the  reference  datums  of  the  tide  gages  were 
not  published. 

Verification.  Mesh  B 

106.  Mesh  B  may  be  viewed  as  an  inset  to  Mesh  A  necessitated  by  two 
factors:  (a)  the  importance  of  riverine/deltaic  systems  to  circulation  in 
portions  of  the  bay  system,  especially  Suisun  Bay,  and  (b)  the  limitation  in 
STUDH  that  allows  only  two-dimensional  elements.  Verification  of  Mesh  B  was 
limited,  since  the  expectation  was  that  by  operation  of  Mesh  B  with  boundary 
conditions  generated  from  Mesh  A  applications  of  RMA-2V,  nearly  identical  re¬ 
sults  should  be  obtained.  This  indeed  was  the  case.  RMA-2V  was  operated  with 
Mesh  A  for  mean  tidal  conditions.  Time-varying  velocity  boundary  conditions 
were  obtained  directly  from  Mesh  A  results  for  the  two-dimensional/one¬ 
dimensional  interfaces  of  the  delta,  Napa  River,  and  Montezuma  and  Suisun 
Sloughs.  Otherwise  Mesh  B  was  operated  with  boundary  conditions  and  model 
coefficients  identical  to  Mesh  A.  Results  of  Mesh  A  and  Mesh  B  were  identical 
based  on  comparison  of  results  at  selected  locations. 

Sediment  Model  Verification 


107.  Concentration  (TSM)  verification  of  the  sediment  model,  STUDH,  was 
performed  to  conditions  and  measurements  during  the  7-8  September  1988  inten¬ 
sive  survey  in  the  Central  Bay  area.  The  verification  data  were  supplemented 
with  limited  surface  TSM  data  obtained  by  the  USGS  during  low  freshwater  in¬ 
flow  periods  during  1969  through  1976  and  1980  (Smith,  Herndon,  and  Harmon 
1979;  Schemel*). 

Model  coefficients  and  parameters 

108.  Several  model  coefficients  and  parameters  required  refinement  dur¬ 
ing  the  verification  process.  These  values  were  selected  to  allow  the  best 
practical  verification  of  natural  sediment  movement.  Special  care  is  neces¬ 
sary  to  deal  with  the  steep  gradients  in  time  and  space  of  natural  sediments 
due  to  resuspension  and  deposition  with  tidal  variations  in  velocities.  At 
some  locations,  temporal  variations  in  TSM  exceeded  one  order  of  magnitude. 

109.  The  diffusion  coefficients  were  specified  at  25  mz/sec  for  the 


*  Personal  communication,  October  1988,  from  L.  Schemel,  US  Geological  Sur¬ 
vey,  Menlo  Park,  CA. 
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x-  and  y-directions .  The  Crank-Nicholson  implicitness  coefficient  9  was  set 
at  0.66,  which  is  the  recommended  value  for  operation  of  STUDH.  (The  Crank- 
Nicholson  scheme  is  used  for  time-stepping  in  STUDH.)  A  value  for  0  of  0.5 
produces  the  most  model  sensitivity  while  a  value  of  0  of  1.0  produces  the 
least  sensitivity.  Increased  values  of  9  smoothed  results,  as  did  increased 
values  of  the  diffusion  coefficients.  The  selected  values  of  0  and  the  dif¬ 
fusion  coefficients  allowed  stable  model  operation  without  excessive  numerical 
and  diffusive  smoothing. 

110.  The  time-step  used  in  RMA-2V  of  37.26  min  (AO  per  24.84-hr  tidal 
cycle)  was  too  large  to  maintain  stable  conditions  in  STUDH.  By  using  a  time- 
step  of  18.63  min  (80  per  tidal  cycle),  STUDH  produced  reasonable  results. 

The  hydrodynamics  required  by  STUDH  at  each  time-step  were  created  through 
postprocessing  of  RMA -2V  results.  In  the  postprocessing,  linear  interpolation 
was  used  to  calculate  the  intermediate  time-step  hydrodynamics  between  RMA-2V 
time-steps . 

111.  Other  model  parameters  and  characteristics  describing  TSM  and 
sediment  beds  are  presented  in  Table  7 .  Boundary  concentrations  were 
estimated  from  the  intensive  field  survey  data  and  historic  surface  TSM  mea¬ 
surements.  TSM  characteristics  and  sediment  bed  characteristics  were  based  on 
previous  laboratory  results  on  San  Francisco  Bay  sediments  (Teeter  1987).  The 
average  settling  velocities  from  separate  settling  tests  conducted  on  the 
September  1988  survey  samples  were  in  the  range  of  0.086  to  0.187  mm/sec. 
Teeter  (1987)  reported  average  settling  velocities  ranging  from  0.437  to 
0.840  mm/sec.  The  process  of  model  verification  gave  satisfactory  results 
with  a  settling  velocity  of  1.0  mm/sec,  which  is  higher  than  the  values  of  the 
September  survey  and  those  reported  in  Teeter  (1987).  However,  -he  size  dis¬ 
tribution  analysis  of  dispersed  samples  from  the  field  survey  indicated  a 
median  particle  size  of  14  /*,  which  is  in  the  silt  range  and  which  substanti¬ 
ates  a  settling  velocity  on  the  order  of  1.0  mm/sec. 

112.  In  the  Gulf  of  the  Faralones,  the  settling  velocity  was  reduced 
beginning  approximately  7  miles  west  of  the  Golden  Gate  Bridge  or  at  the  west¬ 
ern  limits  of  the  Fourfathom  Bank.  The  settling  velocity  was  reduced  from 
0.5  mm/sec  at  the  Fourfathom  Bank  to  0.01  mm/sec  at  the  Farallon  Islands. 

This  was  found  necessary  to  prevent  excessive  deposition  and  loss  of  TSM  con¬ 
centration  in  the  deeper,  slow-moving  gulf  waters.  The  reduced  settling  ve¬ 
locities  in  this  model  area  are  required  because  the  present  numerical  model 
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Table  7 

Cohesive  Sedimentation  Coefficients  and  Model  Parameters 


Coefficient _  Value 


Initial  suspended  sediment 

0. 

.070 

concentration,  kg/m3 

Boundary  concentrations ,  kg/m3 

Ocean 

0, 

.015 

Delta 

0, 

.125 

South  bay 

0. 

.030 

Napa  River 

0 

.060 

Suisun  Slough 

0, 

.060 

Montezuma  Slough 

0 

.040 

Particle  settling  velocity,  mm/sec* 

1, 

.0 

Particle  specific  gravity 

2. 

.65 

Critical  shear  stress  for  deposition  N/m2 

0. 

15 

Critical  shear  stress  for  erosion  N/m2 

0. 

20 

Erosion  rate  constant,  kg/m2/sec 

0. 

0001 

Sediment  bed  initialization 

Noneroding 

Dry  weight  density — newly  deposited  layer,  kg/m3 

450 

Point  source,  Alcatraz  Island,  kg/sec** 

8 

Wind  speed ,  mph 

10 

Wind  fetch  length,  miles 

4 

*  Except  in  ocean  area  (see  text), 

**  Represents  dredged  material  disposal  during  the  period  of  the  September 
1988  study. 
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does  not  incorporate  offshore  currents  known  to  exist  in  the  region,  and  for 
which  inadequate  knowledge  exists  to  quantify  current  magnitudes  accurately. 
Secondly,  at  the  low  concentrations  of  TSM  in  the  gulf,  settling  velocities 
would  be  expected  to  be  reduced.  Generally,  settling  velocities  are  directly 
proportional  to  TSM  concentration,  until  hindered  settling  occurs  at  concen¬ 
trations  much  higher  than  those  of  this  study. 

113.  The  critical  shear  stress  for  erosion  and  deposition  used  in  STUDH 
is  in  the  range  of  values  from  the  University  of  Florida  flume  studies  (Teeter 
1987)  on  San  Francisco  Bay  muds.  However,  in  a  much  earlier  study  (Krone 
1962),  a  minimum  bed  shear  to  allow  deposition  of  0.06  to  0.08  N/m2  was  deter¬ 
mined.  Values  from  Krone  (1962)  used  in  STUDH  produced  poor  results  and  did 
not  provide  adequate  deposition  periods  during  the  tidal  cycle  to  allow  repro¬ 
duction  of  observed  TSM  patterns.  In  STUDH  and  Teeter  (1987),  bed  shear  was 
calculated  by  the  Manning  shear  stress  equation,  whereas  in  Krone  (1962)  the 
smooth-wall  log  velocity  profile  was  used.  Therefore  some  of  the  discrepancy 
in  bed  shear  values  may  be  due  to  the  manner  of  calculation.  (In  ^TUDH,  bed 
shear  can  be  calculated  by  the  smooth-wall  log  velocity  profile  or  the  Manning 
shear  stress  equation.  However,  when  wave-induced  bed  shear  is  important,  as 
it  is  in  San  Francisco  Bay,  a  Jonsson-type  equation  for  wave  and  current 
caused  shear  is  used,  which  incorporates  the  Manning  shear  stress  equation.) 

114.  When  a  Jonsson-type  equation  is  used  to  determine  bed  shear  as  a 
function  of  tidal  currents  and  waves,  a  wind  speed  and  wind  fetch  are  used  to 
calculate  wave  height  and  frequency  in  order  to  determine  wave  orbital  veloc¬ 
ity.  For  the  typical  summer  conditions,  a  daily  average  wind  speed  of  10  mph 
was  selected  with  a  fetch  length  of  4  miles.  Monthly  average  wind  speeds  at 
the  San  Francisco  Airport  for  June  through  October  range  from  9.3  mph  to 
13.9  mph  (NOAA  1985),  which  is  somewhat  greater  than  other  more  inland  sta¬ 
tions.  The  4-mile  fetch  length  is  a  representative  average  for  the  bay  sys¬ 
tem.  This  wind  speed  and  fetch  produced  wave  heights  of  0.5  ft  with  a  period 
of  2  sec.  The  Jonsson-type  equation  has  the  following  general  form: 


1 

\  uom  +  f„  U 

1 

N 

i  U-  +  U 

(U  *  Uom)2 


(9) 
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where 


U*  “  shear  velocity 

fu  =  shear  stress  coefficient  for  waves 
Uom  =  maximum  orbital  velocity  of  waves 
fc  -  shear  stress  coefficients  for  currents 
U  =  vertically  averaged  veloo 

115.  In  the  model  verification,  £.  point  source  was  included  at  the 
Alcatraz  site  (mesh  element  809)  to  approximate  the  level  of  dredged  material 
disposal  that  occurred  during  the  intensive  study  timeframe.  Based  upon  San 
Francisco  District  records,  49,700  cu  yd  of  material  was  disposed  at  the 
Alcatraz  site  in  August  1988  and  44,500  cu  yd  were  disposed  in  September  1988, 
which  are  moderately  low  levels  of  disposal.  All  this  material  came  from 
clamshell  operation.  Using  a  bulk  wet  density  of  1.52  (approximately 

8,000  mg /H)  for  the  in  situ  sediment,  a  15  percent  sand  content,  and  15  per¬ 
cent  material  retention  at  the  Alcatraz  site,  the  average  rate  of  readily 
erodible  cohesive  dredged  material  disposed  at  the  Alcatraz  site  during  August 
and  September  was  8  kg/sec. 

Verification  procedure 

116.  The  STUDH  verification  required  the  following  procedural  steps: 

a.  RMA-2V  was  operated  with  Mesh  A  for  three  tidal  cycles,  not 
including  10  hr  allowed  for  model  spin-up,  with  a  repetitive 
mean  tide  ocean  boundary  condition  and  a  net  delta  outflow  of 
2,500  cfs.  The  average  net  delta  outflow  obtained  from  the 
USGS*  for  the  month  prior  to  the  7-8  September  1988  intensive 
field  survey  was  2,450  cfs  (or  2,500  cfs).  By  the  third  tidal 
cycle,  Mesh  A  hydrodynamics  had  closely  approached  dynamic 
equilibrium  conditions.  (Dynamic  equilibrium  is  the  condition 
in  which  identical  water  levels  and  currents  occur  at  the  same 
tidal  phase  for  successive  tidal  cycles.) 

b.  The  third  tidal  cycle  hydrodynamic  results  from  Mesh  A  pro¬ 
vided  the  necessary  time-varying  velocity  boundary  conditions 
for  Mesh  B. 

c.  RMA-2V  was  operated  with  Mesh  B  for  two  tidal  cycles  with  mean 
tide  conditions.  Since  Mesh  B  did  not  include  the  extensive 
Sacramento-San  Joaquin  Delta  and  used  the  dynamic  velocity 
boundaries  developed  from  Mesh  A,  Mesh  B  approached  dynamic 
equilibrium  more  quickly  than  Mesh  A. 


*  Personal  communication,  November  1988,  from  M.  Taylor,  US  Geological  Sur¬ 
vey,  Sacramento,  CA. 
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d.  The  RMA-2V  hydrodynamics  for  the  second  tidal  cycle  of  opera¬ 
tion  of  Mesh  B  were  sent  through  appropriate  postprocessing 
codes  to  develop  the  hydrodynamics  necessary  to  operate  STUDH. 

e.  Beginning  with  constant  TSM  value  of  0.07  kg/m3  and  the  other 
coefficients  shown  in  Table  7,  STUDH  was  operated  for  three 
repetitive  tidal  cycles  with  the  hydrodynamics  developed  for 
Mesh  B.  The  three  tidal  cycles  allowed  ample  time  for  the 
solution  to  adjust  from  initial  conditions  and  to  develop 
depositional  beds,  erosion/deposition  patterns,  and  TSM  pat¬ 
terns  both  spatially  and  temporally. 

f.  Using  the  ending  hydrodynamics  of  the  third  tidal  cycle  from 
step  a,  RMA-2V  was  operated  using  Mesh  A  for  the  verification 
period  tidal  cycle  of  8  a.m.,  7  September,  to  9  a.m. , 

8  September  1988.  (The  mean  tide  of  steps  a  through  e  was 
oriented  in  time  so  that  it  could  be  directly  phased  into  the 
verification  period  tide  in  a  smooth,  consistent  manner.) 

g.  The  Mesh  A  results  for  the  7-8  September  tide  were  used  to 
generate  time-varying  velocity  boundary  conditions  for  Mesh  B. 

h.  RMA-2V  was  operated  with  Mesh  B  for  the  7-8  September  tide. 

i. .  The  hydrodynamic  results  from  Mesh  B  were  sent  through  post¬ 

processing  to  develop  the  hydrodynamics  to  operate  STUDH  for 
7-8  September. 

j. .  Using  the  ending  TSM  concentrations  and  sediment  bed  structure 

from  the  third  mean  tide  cycle,  STUDH  was  operated  for  the 
verification  period  of  7-8  September  1988. 

Verification  results 

117.  The  results  from  RMA-2V  for  Mesh  B  simulation  of  the  verification 
period  of  8  a.m.,  7  September,  to  9  a.m.,  8  September  (model  hour  50  to  75), 
are  displayed  with  field-measured  velocities  in  Plates  12-14.  (Figure  3  dis¬ 
plays  the  field  sample  sites  used  in  verification.)  The  field  measurements 
were  weight-averaged  vertically  (see  next  paragraph)  to  obtain  values  for 
comparison  with  RMA-2V  results.  Good  comparisons  of  RMA-2V  and  field  measure¬ 
ment  were  obtained  at  all  locations  with  the  exceptions  of  the  very  shallow 
water,  sta  4D  (depth)  and  to  a  lesser  extent  sta  4C  (depth).  The  exact  cause 
of  the  poorer  comparisons  at  these  shallow  stations  is  unknown,  but  some  dis¬ 
crepancy  is  probably  the  result  of  high  wind  and  wave  action  corrupting  field 
measurements  at  these  shallow,  lower  velocity  stations.  The  phasing  of  numer¬ 
ical  results  is  in  good  agreement  with  the  physical  data  at  til  1  Sta  tions ,  and 
velocity  magnitude  comparisons  are  good  especially  at  sta  3A,  3B,  and  3C. 

118.  The  results  from  STUDH  for  the  verification  period  8  a.m., 

7  September,  through  9  a.m.,  8  September,  are  displayed  with  the  field- 
measured  verification  data  in  Plates  15-21.  (Figure  3  displays  the  field 
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sample  sites  used  in  verification.)  The  verification  data  from  the  field 
measurements  were  vertically  averaged  for  comparison  with  the  vertically 
averaged  STUDH  results.  The  field  data  averaging  process  involved  weighting 
to  approximate  the  percent  of  the  vertical  water  column  represented  by  a 
value.  For  stations  with  three  vertical  measurements,  the  vertically  averaged 
value  X  was  determined  as 

X  =  0.25  XB  +  0.50  Xm  +  0.25  Xb  (10) 

where 

Xfl  **  near-surface  value 
X„,  -  middepth  value 
XB  “  near-bottom  value 

For  the  situation  with  five  vertical  measurements,  then 


X  =  0.1Xa  +  0 . 25Xx  +  0.3\  +  0 . 25X3  +  0.1X,, 


(11) 


where 

Xx  =  one-fourth  depth  value 

X3  “  three-fourths  depth  value 

A  simple  arithmetic  mean  was  not  used  in  Equations  10  and  11  because  the  near¬ 
surface  and  near-bottom  measurements  are  not  representative  of  the  same  per¬ 
centage  of  the  water  column  as  an  interior  measurement  point.  The  numerical 
results  are  generally  satisfactory  considering  the  following:  potential 
biases  and  errors  induced  from  sampling  and  the  high  variability  of  natural 
TSM  values  in  San  Francisco  Bay  in  space  and  time  resulting  from  the  interac¬ 
tion  of  periods  of  deposition  followed  by  periods  of  resuspension.  In  the 
numerical  and  prototype  TSM  values,  the  fluctuations  within  tidal  cycles  were 
pronounced.  Generally  during  slack  water,  TSM  values  decreased  only  to  in¬ 
crease  rapidly  in  response  to  velocity  increase  (resulting  in  a  resuspension 
of  material) .  The  numerical  model  results  generally  mimicked  the  measured 
patterns  and  responded  similar  to  current  speed  (and  thus  bed  shear  and  sedi¬ 
mentation  (deposition  and  resuspension))  variations.  A  notable  exception  was 
the  numerical  model's  inability  to  reproduce  the  pronounced  magnitude  of  in¬ 
crease  in  TSM  concentration  occurring  during  the  strong  ebb  current  from  hhw 
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to  llw  at  approximately  hour  18  on  the  plots  or  approximately  2  a.m.  , 

8  September;  for  example,  see  sta  1C,  3A,  and  3B.  This  temporary  surge  in 
concentration  was  adequately  reproduced  at  shallow  sta  4A  and  4D.  The  com¬ 
plexity  of  sediment  movement  in  the  bay  system  is  indicated  by  survey  measure¬ 
ments  showing  no  significant  surge  during  the  strong  ebb  current  at  sta  3C, 
though  a  sharp  concentration  rise  was  measured  at  approximately  11  p.m., 

7  September,  near  slack  water,  while  adjacent  sta  3A  and  3B  showed  the  ebb 
current  response.  A  general  weakness  of  the  numerical  results  is  in  reproduc¬ 
ing  the  magnitude  of  the  response  of  the  prototype  system  to  variations  in  bed 
shear  manifested  through  the  temporal  changes  (acceleration/deceleration)  in 
tidal  currents.  The  numerical  simulation  also  had  some  difficulty  in  repro¬ 
ducing  the  lateral  variation  in  TSM  across  Range  1,  sta  1A,  IB,  and  1C.  Rela¬ 
tively  good  simulations  were  obtained  at  sta  IB,  2A,  2B,  2C,  4A,  and  4B. 

119.  As  a  second  less  exacting  verification,  historical  USGS  TSM  data 
were  used  (Smith,  Herndon,  and  Harmon  1979;  Schemel*) .  These  data  were  taken 
along  the  deep  ship  channels  transecting  the  bay  system  from  the  delta  to  the 
end  of  South  Bay  and  were  generally  taken  at  a  depth  of  2  ft  (surface  read¬ 
ings).  Based  on  the  WES  intensive  survey  data,  an  average  factor  by  which  to 
multiply  a  surface  concentration  to  estimate  the  vertical  average  value  is 
1.6.  This  factor  was  determined  by  dividing  the  average  of  all  surface  TSM 
concentrations  for  the  7-8  September  intensive  survey  into  the  average  verti¬ 
cally  averaged  value  for  each  deepwater  station  (excluding  shallow  sta  4C  and 
4D) .  The  factor  ranged  from  1.2  to  2.0.  While  it  is  not  exacting  to  compare 
numerically  determined  TSM  concentrations  from  one  tide  condition  and  one  net 
delta  outflow  to  field  measurements  taken  during  various  tide  conditions  from 
spring  to  neap  and  for  a  wide  range  of  low  to  extremely  low  net  delta  out¬ 
flows,  the  exercise  did  indicate  whether  the  model  replicated  the  generally 
observed  trends  in  the  natural  system  (Table  8).  The  stations  for  comparison 
had  a  minimum  of  17  surface  TSM  measurements  taken  during  the  summers  of  1972, 
1973,  1976,  and  summer  and  fall  of  1980,  which  were  periods  with  low  monthly 
average  delta  outflows  ranging  from  3,000  to  15,000  cfs.  The  numerical  model 
minimum,  maximum,  and  mean  values  for  nodes  close  to  the  USGS  stations  gener¬ 
ally  compared  very  favorably  with  the  field  measurements.  In  Suisun  Bay 
(sta  6),  the  numerical  results  appeared  to  be  somewhat  low,  whereas  in  South 


*  Op.  cit. 
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Bay  (sta  24),  the  numerical  results  were  too  high.  For  the  remaining  stations 
in  Central  and  San  Pablo  Bay,  the  numerical  and  field  results  compared 
favorably . 
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PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


120.  This  study  was  limited  in  that  only  a  late  summer/early  fall  low 
freshwater  inflow  condition  was  considered.  The  subsequent  study  phase  will 
include  the  winter/spring  high  freshwater  inflow  period  so  that  seasonal 
variations  within  the  system  can  be  determined.  Upon  the  completion  of  that 
phase  of  this  study,  the  results  of  the  two  phases  will  be  integrated  and  an 
explanation  of  the  fate  of  dredged  material  will  be  presented. 

Field  Data  Conclusions 


121.  The  field  data  collected  during  the  low  freshwater  inflow  period 
of  September  1988  indicated  that  Central  Bay  had  near-oceanic  salinities,  was 
vertically  well  mixed,  and  had  moderately  low  total  suspended  material  (TSM) . 
Net  delta  outflow  of  fresh  water  into  the  bay  system  was  only  about  2,500  cfs 
during  this  period.  TSM  at  the  Golden  Gate  survey  range  averaged  about 

16  mg/i,  increased  upstream  to  34  mg/i  at  the  entrance  to  South  Bay,  and  in¬ 
creased  to  45  mg /Si  at  the  Richmond  Bridge  near  the  entrance  to  San  Pablo  Bay. 
Peak  TSM  values  were  greater  than  average  values  by  a  factor  of  about  three. 

122.  Tidal  fluxes  of  TSM  were  large,  commensurate  with  the  tidal  prism 
of  this  system.  Observed  tidal  fluxes  were  roughly  30,000  metric  tons  per 
lunar  day  at  the  Golden  Gate  range,  19,000  metric  tons  per  lunar  day  at  the 
entrance  to  the  South  Bay,  and  44,000  metric  tons  per  lunar  day  at  Richmond. 

123.  Despite  the  genial  tidal  characteristic  of  Central  Bay,  which  fa¬ 
vors  ebb-dominated  flows,  net  landward  TSM  fluxes  were  observed  through  the 
ranges  at  the  Golden  Gate,  and  the  entrance  to  South  Bay.  These  landward 
fluxes  had  relatively  small  magnitudes,  5-10  percent  of  the  tidal  fluxes.  The 
station  data  suggest  possible  clockwise  net  TSM  flux  for  that  part  of  Central 
Bay  where  estimates  could  be  made.  The  net  fluxes  were  south  along  the  east 
side  of  the  bay  near  the  Richmond  Bridge,  Berkeley  shallows,  and  into  South 
Bay.  The  net  flux  was  northward  on  the  west  side  of  the  entrance  to  South 
Bay. 

124.  The  overall  flux  balance  for  Central  Bay  was  not  determined.  The 
flux  at  sta  1C  at  the  Richmond  Bridge  was  seaward,  but  a  large  portion  of  the 
range  cross  section  was  not  sampled. 

125.  The  phenomenon  responsible  for  the  largest  fraction  of  the  TSM 
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flux  was  fluctuation  in  TSM  during  high  current  speeds,  called  tidal  pumping. 
Such  fluctuations  were  most  pronounced  near  the  bottom.  TSM  fluctuations  were 
apparently  more  rapid  than  the  1-hr  sampling  rate,  and  were  therefore  not  well 
resolved. 

126.  The  question  remains  how  well  the  intensive  survey  was  representa¬ 
tive  of  low  freshwater  inflow  conditions .  Tidal  and  wind  conditions  were 
about  normal.  There  were  no  strong  variations  in  salinity  or  TSM  conditions 
during  the  2-week  period  surrounding  the  intensive  survey.  However,  the  vari¬ 
ability  in  fluxes  over  days  or  weeks  is  not  known,  and  therefore  how  well  the 
intensive  survey  represents  seasonal  conditions  remains  an  open  question. 

127.  Data  from  previous  published  studies  were  used  in  an  attempt  to 
determine  near-bed  tidal  dominance.  Although  results  were  not  always  con¬ 
sistent,  a  complex  pattern  with  both  ebb-  and  flood-dominant  areas  across  most 
cross  sections  emerged.  The  area  near  the  Alcatraz  disposal  site  was 
indicated  to  be  ebb  dominated. 

Dye  Test  Conclusions 

128.  Previous  physical  model  dye  test  data  from  the  San  Francisco  Bay- 
Delta  Physical-Hydraulic  Model  (SFBM)  were  reanalyzed  to  determine  dispersion 
of  dissolved  materials.  Although  dyes  in  general  are  a  poor  analog  for  sus¬ 
pended  sediments ,  they  are  representative  of  nonsettling  wash  load  and  are 
affected  by  the  same  flow  and  mixing  regimes  that  influence  sediment  trans¬ 
port.  Flushing  rates  and  residence  times  were  estimated  for  five  segments  in 
the  Central  Bay  area.  Results  indicated  that  most  segments  had  hydraulic 
residence  times  of  1-2  lunar  days,  except  for  shallow  areas,  which  had  resi¬ 
dence  times  of  4-7  lunar  days.  A  maximum  concentration  increase  of  about 

2  mg/X  for  any  segment  was  estimated  for  a  continuous  discharge  of 
10,000  metric  tons  per  tidal  cycle  near  the  Alcatraz  disposal  site. 

Numerical  Modeling  Conclusions 

129.  The  model  RMA-2V,  which  predicts  two-dimensional  vertically  aver¬ 
aged  velocities,  was  successfully  verified  for  the  San  Francisco  Bay  system  to 
water  levels  developed  from  harmonic  analysis  of  prototype  tide  data  and  to 
physical  model  velocities  and  water  levels  from  the  SFBM. 
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130.  The  sediment  transport  model,  STUDH,  which  uses  the  hydrodynamic 
results  from  RMA-2V,  was  verified  to  field-measured  total  suspended  matter 
(TSM)  concentrations  from  an  intensive  survey  conducted  for  25  hr  during 
7-8  September  1988  and  to  historical  field  data  gathered  by  the  US  Geological 
Survey  (USGS).  STUDH  was  only  moderately  successful  in  replicating  the  highly 
variable  TSM  values  obtained  during  the  intensive  survey.  However,  the  model 
results  did  compare  favorably  with  the  less  exacting  ranges  of  the  historical 
USGS  data.  Generally,  the  model  results  did  not  completely  reproduce  the 
short  pulses  of  high  TSM  concentrations  of  the  field  survey,  which  frequently, 
though  not  always,  coincided  with  periods  of  high  tidal  velocities  and  bed 
shear. 


Recommendations 


131.  Several  areas  of  field  data  paucity  were  identified  during  this 
study.  Further  coordinated  field  investigations  are  recommended  to  improve 
the  general  understanding  and  predictive/diagnostic  modeling  of  suspended 
sediment  transport  in  Central  Bay. 

a.  Map  surficial  bottom  sediment  size,  metals  content,  and 
physico-chemical  characteristics  related  to  fine-sediment 
erodibility  seasonally  over  an  annual  cycle,  and,  if  possible, 
over  storm  events. 

b.  Analyze  historical  surveys  to  identify  long-term  deposition 
and  erosion  trends. 

c.  Deploy  near-bed  instrument  arrays  at  various  positions  in 
Central  Bay  to  obtain  high  temporal  descriptions  of  suspended 
sediment  pulses,  and  associated  flow  conditions.  Simulta¬ 
neously,  investigate  local  suspended  sediments  at  high  spatial 
resolution  using  towed  sensors. 

132.  To  complete  the  analysis  of  long-term  transport,  additional  data 
sets  are  needed  to  quantify  transport  conditions  and  for  use  in  future  sedi¬ 
ment  model  verification.  Such  data  sets  should  include,  at  the  minimum,  a 
high-flow  wet-weather  survey,  and  possibly  a  winter  frontal  passage  condition. 
Duplicate  data  sets  would  be  desirable.  Areas  of  long-term  deposition  and 
erosion  within  Central,  San  Pablo,  and  South  Bays  should  be  identified  and 
quantified  by  comparisons  of  historical  bathymetric  surveys,  radioisotope 
tracer  studies,  and/or  grain  size  characterizations. 

133.  Several  areas  of  potential  numerical  model  modification, 


76 


improvement,  or  redirection  are  recommended  for  future  sediment  modeling  of 
the  Alcatraz  site  and  San  Francisco  Bay: 

a.  The  ocean  boundary  in  the  hydrodynamic  model,  RMA-2V,  should 
be  modified  to  produce  an  appropriate,  seasonal  longshore 
current  in  the  Gulf  of  the  Farallones.  This  will  more  accu¬ 
rately  define  the  exchange  of  waters  between  the  gulf  and  San 
Francisco  Bay.  The  nearshore  currents  off  the  California 
coast  can  be  very  complex,  seasonally  varying,  and  subject  to 
anomalies  most  often  attributed  to  eddies,  meanders,  and 
counterflows.  The  California  current  is  a  slow  southward 
current  generally  close  to  shore  during  late  winter  through 
early  fall.  A  northward  countercurrent,  the  Davidson  Current, 
persists  nearshore  generally  in  late  fall  and  early  winter. 
There  is  some  limited  evidence  of  a  nearshore  poleward  flow 
throughout  the  year  along  the  central  California  coast  that  is 
not  resolved  by  normal  hydrographic  data,  which  is  taken  too 
far  offshore. 

b.  The  present  numerical  mesh  should  be  extended  to  include  the 
southern  extremity  of  South  Bay,  below  Dumbarton  Bridge,  and 
to  refine  the  resolution  in  the  immediate  vicinity  of  the 
Alcatraz  site.  The  mesh  extension  will  directly  resolve  a 
difficult  boundary  condition  and  the  mesh  refinement  will 
facilitate  the  accuracy  and  stability  of  the  material  injec¬ 
tion  at  the  disposal  site. 

c.  Settling  velocity  should  be  included  as  a  function  of  TSM  con¬ 
centration  in  STUDH.  Such  a  relationship  has  been  noted  in 
settling  tests  on  sediments  from  San  Francisco  Bay,  as  well  as 
other  estuarine  systems. 

d.  Wind-wave  induced  bed  shear  should  be  incorporated  more  accu¬ 
rately  as  a  function  of  wind  direction  and  actual  fetch  in 
STUDH.  Because  wind  direction  and  speed  varies  both  tempor¬ 
ally  and  spatially  over  the  San  Francisco  Bay  system,  this  may 
be  an  important  consideration  and  will  require  the  application 
of  a  wind  model. 

e.  The  wind-surface  stress  and  density  gradient  options  in  opera¬ 
tion  of  RMA-2V  should  be  included,  especially  to  determine  the 
influence  on  STUDH  verification  results  to  the  intensive 
survey  data. 

f .  Improved  boundary  conditions  should  be  obtained,  especially 
TSM  values,  for  the  Gulf  of  the  Faralones.  This  can  be 
obtained  by  a  more  thorough  literature  search  and/or  field 
measurements . 

g.  The  benefit  of  having  multiple  grain  sizes,  essentially  silt 
and  <-lay,  for  representation  of  fine-grained  material  in  STUDH 
should  be  investigated.  The  observed  rapid  temporal  fluctua¬ 
tions  in  field  TSM  may  be  the  result  of  a  silt  fraction  that 
settles  rapidly.  Because  of  the  computational  requirements, 
this  option  becomes  viable  only  when  a  large  mainframe  or 
supercomputer  is  available.  The  WES  supercomputer  became 
available  during  January  1990. 
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134.  The  viability  of  three-dimensional  hydrodynamic  and  sediment 
transport  numerical  modeling  should  be  investigated.  A  coarse  mesh  for 
limited  three-dimensional  hydrodynamic-salinity  simulations  was  used  on  the 
John  F.  Baldwin  Phase  III  Ship  Channel  Project,  and  the  mesh  is  available  for 
adaptation  as  a  three-dimensional  sediment  mesh.  However,  three-dimensiom.1 
sediment  modeling  is  not  presently  an  off-the-shelf  application,  and  must  be 
approached  with  much  forethought. 
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APPENDIX  B:  THE  TABS -2  SYSTEM 


1.  TABS -2  is  a  collection  of  generalized  computer  programs  and  utility 
codes  integrated  into  a  numerical  modeling  system  for  studying  two-dimensional 
hydrodynamics,  sedimentation,  and  transport  problems  in  rivers,  reservoirs, 
bays,  and  estuaries.  A  schematic  representation  of  the  system  is  shown  in 
Figure  Bl.  It  can  be  used  either  as  a  stand-alone  solution  technique  or  as  a 
step  in  the  hybrid  modeling  approach.  The  basic  concept  is  to  calculate 
water-surface  elevations,  current  patterns,  sediment  erosion,  transport  and 
deposition,  the  resulting  bed  surface  elevations,  and  the  feedback  to  hydrau¬ 
lics.  Existing  and  proposed  geometry  can  be  analyzed  to  determine  the  impact 
on  sedimentation  of  project  designs  and  to  determine  the  impact  of  project 
designs  on  salinity  and  on  the  stream  system.  The  system  is  described  in  de¬ 
tail  by  Thomas  and  McAnally  (1985) . 

2.  The  three  basic  components  of  the  system  are  as  follows: 

a.  "A  Two-Dimensional  Model  for  Free  Surface  Flows,"  RMA-2V. 

b.  "Sediment  Transport  in  Unsteady  2-Dimensional  Flows,  Horizontal 
Plane,"  STUDH. 

c.  "Two-Dimensional  Finite  Element  Program  for  Water  Quality," 
RMA-4. 

3.  RMA-2V  is  a  finite  element  solution  of  the  Reynolds  form  of  the 
Navier-Stokes  equations  for  turbulent  flows.  Friction  is  calculated  with 
Manning's  equation  and  eddy  viscosity  coefficients  are  used  to  define  the 
turbulent  losses.  A  velocity  form  of  the  basic  equation  is  used  with  side 
boundaries  treated  as  either  slip  or  static.  The  model  automatically  recog¬ 
nizes  dry  elements  and  corrects  the  mesh  accordingly.  Boundary  conditions  may 
be  water-surface  elevations,  velocities,  or  discharges  and  may  occur  inside 
the  mesh  as  well  as  along  the  edges. 


TABS-2 


Figure  Bl.  TABS-2  schematic 
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4.  The  sedimentation  model,  STUDH,  solves  the  convection- diffusion 
equation  with  bed  source  terms.  These  terms  are  structured  for  either  sand  or 
cohesive  sediments.  The  Ackers-White  (1973)  procedure  is  used  to  calculate  a 
sediment  transport  potential  for  the  sands  from  which  the  actual  transport  is 
calculated  based  on  availability.  Clay  erosion  is  based  on  work  by  Parthen- 
iades  (1962)  and  Ariathurai  and  the  deposition  of  clay  utilizes  Krone's  equa¬ 
tions  (Ariathurai,  MacArthur,  and  Krone  1977).  Deposited  material  forms 
layers,  as  shown  in  Figure  B2,  and  bookkeeping  allows  up  to  10  layers  at  each 
node  for  maintaining  separate  material  types,  deposit  thickness,  and  age.  The 
code  uses  the  same  mesh  as  RMA-2V. 

5.  Salinity  calculations,  RMA-4,  are  made  with  a  form  of  the 
convective -diffusion  equation  which  has  general  source -sink  terms.  Up  to 
seven  conservative  substances  or  substances  requiring  a  decay  term  can  be 
routed.  The  code  uses  the  same  mesh  as  RMA-2V. 

6.  Each  of  these  generalized  computer  codes  can  be  used  as  a  stand¬ 
alone  program,  but  to  facilitate  the  preparation  of  input  data  and  to  aid  in 
analyzing  results,  a  family  of  utility  programs  was  developed  for  the  follow- 


ing  purposes : 

a . 

Digitizing 

b. 

Mesh  generation 

c. 

Spatial  data  management 

d. 

Graphical  output 

e . 

Output  analysis 

f. 

File  management 

&■ 

Interfaces 

h. 

Job  control  language 

Finite  Element  Modeling 

7.  The  TABS-2  numerical  models  used  in  this  effort  employ  the  finite 
element  method  to  solve  the  governing  equations.  To  help  those  who  are  un¬ 
familiar  with  the  method  to  better  understand  this  report,  a  brief  description 
of  the  method  is  given  here. 

8.  The  finite  element  method  approximates  a  solution  to  equations  by 
dividing  the  area  of  interest  into  smaller  subareas,  which  are  called  ele¬ 
ments.  The  dependent  variables  (e.g.,  water-surface  elevations  and  sediment 
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concentrations)  are  approximated  over  each  element  by  continuous  functions 
which  interpolate  in  terms  of  unknown  point  (node)  values  of  the  variables. 

An  error,  defined  as  the  deviation  of  the  approximation  solution  from  the  cor¬ 
rect  solution,  is  minimized.  Then,  when  boundary  conditions  are  imposed,  a 
set  of  solvable  simultaneous  equations  is  created.  The  solution  is  continuous 
over  the  area  of  interest. 

9.  In  one -dimensional  problems,  elements  are  line  segments.  In  two- 
dimensional  problems,  the  elements  are  polygons,  usually  either  triangles  or 
quadrilaterals.  Nodes  are  located  on  the  edges  of  elements  and  occasionally 
inside  the  elements.  The  interpolating  functions  may  be  linear  or  higher 
order  polynomials.  Figure  B2  illustrates  a  quadrilateral  element  with  eight 
nodes  and  a  linear  solution  surface  where  F  is  the  interpolating  function. 

10.  Most  water  resource  applications  of  the  finite  element  method  use 
the  Galerkin  method  of  weighted  residuals  to  minimize  error.  In  this  method 
the  residual,  the  total  error  between  the  approximate  ai.J  correct  solutions, 
is  weighted  by  a  function  that  is  identical  with  the  interpolating  function 
and  then  minimized.  Minimization  results  in  a  set  of  simultaneous  equations 
in  terms  of  nodal  values  of  the  dependent  variable  (e.g.  water-surface  eleva¬ 
tions  or  sediment  concentration).  The  time  portion  of  time -dependent  problems 
can  be  solved  by  the  finite  element  method,  but  it  is  generally  more  efficient 
to  express  derivatives  with  respect  to  time  in  finite  difference  form. 

The  Hydrodynamic  Model.  RMA-2V 


Applications 

11.  This  program  is  designed  for  far-field  problems  in  which  vertical 
accelerations  are  negligible  and  the  velocity  vectors  at  a  node  generally 
point  in  the  same  directions  over  the  entire  depth  of  the  water  column  at  any 
instant  of  time.  It  expects  a  homogeneous  fluid  with  a  free  surface.  Both 
steady  and  unsteady  state  problems  can  be  analyzed.  A  surface  wind  stress  can 
be  imposed. 

12.  The  program  has  been  applied  to  calculate  flow  distribution  around 
islands;  flow  at  bridges  having  one  or  more  relief  openings,  in  contracting 
and  expanding  reaches,  into  and  out  of  off -channel  hydropower  plants,  at  river 
junctions,  and  into  and  out  of  pumping  plant  channels;  and  general  flow  pat¬ 
terns  in  rivers,  reservoirs,  and  estuaries. 
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Limitations 


13.  This  program  is  not  designed  for  near-field  problems  where  flow- 
structure  interactions  (such  as  vortices,  vibrations,  or  vertical  accelera¬ 
tions)  are  of  interest.  Areas  of  vertically  stratified  flow  are  beyond  this 
program's  capability  unless  it  is  used  in  a  hybrid  modeling  approach.  It  is 
two-dimensional  in  the  horizontal  plane,  and  zones  where  the  bottom  current  is 
in  a  different  direction  from  the  surface  current  must  be  analyzed  with  con¬ 
siderable  subjective  judgment  regarding  long-term  energy  considerations.  It 
is  a  free-surface  calculation  for  subcritical  flow  problems. 

Governing  equations 

14.  The  generalized  computer  program  RMA-2V  solves  the  depth- integrated 
equations  of  fluid  mass  and  momentum  conservation  in  two  horizontal  direc¬ 
tions  .  The  form  of  the  solved  equations  is 
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where 

h  -  depth 

u,v  =  velocities  in  the  Cartesian  directions 
x,y,t  =  Cartesian  coordinates  and  time 
p  =  density 
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£  -  eddy  viscosity  coefficient,  for  xx  -  normal  direction  on 

x-axis  surface;  yy  -  normal  direction  on  y-axis  surface;  xy 
and  yx  =  shear  direction  on  each  surface 

g  -  acceleration  due  to  gravity 

a  -  elevation  of  bottom 

n  -  Manning's  n  value 

1.486  -  conversion  from  SI  (metric)  to  non-SI  units 
f  -  empirical  wind  shear  coefficient 
-  wind  speed 
rj>  -  wind  direction 
w  “  rate  of  earth's  angular  rotation 
<f>  -  local  latitude 

15.  Equations  Bl,  B2,  and  B3  are  solved  by  the  finite  element  method 
using  Galerkin  weighted  residuals.  The  elements  may  be  either  quadrilaterals 
or  triangles  and  may  have  curved  (parabolic)  sides.  The  shape  functions  are 
quadratic  for  flow  and  linear  for  depth.  Integration  in  space  is  performed  by 
Gaussian  integration.  Derivatives  in  time  are  replaced  by  a  nonlinear  finite 
difference  approximation.  Variables  are  assumed  to  vary  over  each  time  inter¬ 
val  in  the  form 


/( t)  »  /( 0)  +  at  +  btC  tp  <  t  <  t  (B4) 

which  is  differentiated  with  respect  to  time,  and  cast  in  finite  difference 
form.  Letters  a  ,  b  ,  and  c  are  constants.  It  has  been  found  by  experi¬ 
ment  that  the  best  value  for  c  is  1.5  (Norton  and  King  1977). 

16.  The  solution  is  fully  implicit  and  the  set  of  simultaneous  equa¬ 
tions  is  solved  by  Newton-Raphson  iteration.  The  computer  code  executes  the 
solution  by  means  of  a  front- type  solver  that  assembles  a  portion  of  the 
matrix  and  solves  it  before  assembling  the  next  portion  of  the  matrix.  The 
front  solver's  efficiency  is  largely  independent  of  bandwidth  and  thus  does 
not  require  as  much  care  in  formation  of  the  computational  mesh  as  do  tradi¬ 
tional  solvers. 

17.  The  code  RMA-2V  is  based  on  the  earlier  version  RMA-2  (Norton  and 
King  1977)  but  differs  from  it  in  several  ways.  It  is  formulated  in  terms  of 
velocity  (v)  instead  of  unit  discharge  (vh) ,  which  improves  some  aspects  of 
the  code's  behavior;  it  permits  drying  and  wetting  of  areas  within  the  grid; 
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and  it  permits  specification  of  turbulent  exchange  coefficients  in  directions 
other  than  along  the  x-  and  z-axes.  Tor  a  more  complete  description,  see 
Appendix  F  of  Thomas  and  McAnally  (1985) . 

The  Sediment  Transport  Model.  STUDH 


Applications 

18.  STUDH  can  be  applied  to  clay  and/or  sand  bed  sediments  where  flow 
velocities  can  be  considered  t^o- dimensional  (i.e.,  the  speed  and  direction 
can  be  satisfactorily  represented  as  a  depth- averaged  velocity).  It  is  useful 
for  both  deposition  and  erosion  studies  and,  to  a  limited  extent,  for  stream 
width  studies.  The  program  treats  two  categories  of  sediment:  noncohesive, 
which  is  referred  to  as  sand  here,  and  cohesive,  which  is  referred  to  as  clay. 
Limitations 

19.  Both  clay  and  sand  may  be  analyzed,  but  the  model  considers  a 
single,  effective  grain  size  for  each  and  treats  each  separately.  Fall  veloc¬ 
ity  must  be  prescribed  along  with  the  water-surface  elevations,  x-velocity, 
y-velocity,  diffusion  coefficients,  bed  density,  critical  shear  stresses  for 
erosion,  erosion  rate  constants,  and  critical  shear  stress  for  deposition. 

20.  Many  applications  cannot  use  long  simulation  periods  because  of 
their  computation  cost.  Study  areas  should  be  made  as  small  as  possible  to 
avoid  an  excessive  number  of  elements  when  dynamic  runs  are  contemplated  yet 
must  be  large  enough  to  permit  proper  posing  of  boundary  conditions.  The  same 
computation  time  interval  must  be  satisfactory  for  both  the  transverse  and 
longitudinal  flow  directions. 

21.  The  program  does  not  compute  water-surface  elevations  or  veloci¬ 
ties;  therefore  these  data  must  be  provided.  For  complicated  geometries,  the 
numerical  model  for  hydrodynamic  computations,  RMA-2V,  is  used. 

Governing  equations 

22.  The  generalized  computer  program  STUDH  solves  the  depth- integrated 
convection-dispersion  equation  in  two  horizontal  dimensions  for  a  single  sedi¬ 
ment  constituent.  For  a  more  complete  description,  see  Appendix  G  of  Thomas 
and  McAnally  (1985) .  The  form  of  the  solved  equation  is 
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where 


C  -  concentration  of  sediment 
u  -  depth- integrated  velocity  in  x-direction 
v  -  depth- integrated  velocity  in  y-direction 
-  dispersion  coefficient  in  x-direction 

Dy  -  dispersion  coefficient  in  y-direction 

a ^  -  coefficient  of  concentration- dependent  source/sink  term 
“  coefficient  of  source/sink  term 

23.  The  source/sink  terms  in  Equation  B5  are  computed  in  routines  that 
treat  the  interaction  of  the  flow  and  the  bed.  Separate  sections  of  the  code 
handle  computations  for  clay  bed  and  sand  bed  problems. 

Sand  transport 

24.  The  source/sink  terms  are  evaluated  by  first  computing  a  potential 
sand  transport  capacity  for  the  specified  flow  conditions,  comparing  that 
capacity  with  the  amount  of  sand  actually  being  transported,  and  then  eroding 
from  or  depositing  to  the  bed  at  a  rate  that  would  approach  the  equilibrium 
value  after  sufficient  elapsed  time. 

25.  The  potential  sand  transport  capacity  in  the  model  is  computed  by 
the  method  of  Ackers  and  White  (1973) ,  which  uses  a  transport  power  (work 
rate)  approach.  It  has  been  shown  to  provide  superior  results  for  transport 
under  steady- flow  condi^ons  (White,  Milli,  and  Crabbe  1975)  and  for  combined 
waves  and  currents  (Swart  1976) .  Flume  tests  at  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  have  shown  that  the  concept  is  valid  for  transport  by 
estuarine  currents. 

26.  The  total  load  transport  function  of  Ackers  and  White  is  based  upon 
a  dimensionless  grain  size 
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where 

D  =  sediment  particle  diameter 
s  =  specific  gravity  of  the  sediment 
v  =  kinematic  viscosity  of  the  fluid 
and  a  sediment  mobility  parameter 
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where 

t  -  total  boundary  shear  stress 

n'  -  a  coefficient  expressing  the  relative  importance  of  bed-load  and 
suspended- load  transport,  given  in  Equation  A9 

t'  **  boundary  surface  shear  stress 

The  surface  shear  stress  is  that  part  of  the  total  shear  stress  which  is  due 
to  the  rough  surface  of  the  bed  only,  i.e.,  not  including  that  part  due  to  bed 
forms  and  geometry.  It  therefore  corresponds  to  that  shear  stress  that  the 
flow  would  exert  on  a  plane  bed. 

27.  The  total  sediment  transport  is  expressed  as  an  effective 
concentration 
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where  U  is  the  average  flow  speed,  and  for  1  <  D  <60 
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For  D  <60 
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A  -  0.17 
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C  -  0.025  (B15) 

m  -  1.5  (Bib) 


28.  Equations  B6-B16  result  in  a  potential  sediment  concentration  . 
This  value  is  the  depth -averaged  concentration  of  sediment  that  will  occur  if 
an  equilibrium  transport  rate  is  reached  with  a  nonlimited  supply  of  sediment. 
The  rate  cf  sediment  deposition  (or  erosion)  is  then  computed  as 
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where 

C  =  present  sediment  concentration 

t  -  time  constant 
c 

For  deposition,  the  time  constant  is 


larger  of 


At 

or 

Cdh 


and  for  erosion  it  is 


At 


tc  =  larger  of 
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where 

At  -  computational  time -step 

=  response  time  coefficient  for  deposition 
V  ”  sediment  settling  velocity 
Ce  =  response  time  coefficient  for  erosion 
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The  sand  bed  has  a  specified  initial  thickness  which  limits  the  amount  of  ero¬ 
sion  to  that  thickness. 

Cohesive  sediments  transport 

29.  Cohesive  sediments  (usually  clays  and  some  silts)  are  considered  to 
be  depositional  if  the  bed  shear  stress  exerted  by  the  flow  is  less  than  a 
critical  value  .  When  that  value  occurs,  the  deposition  rate  is  given  by 
Krone's  (1962)  equation 
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where 


S  -  source  term 

Vg  -  fall  velocity  of  a  sediment  particle 
h  -  flow  depth 

C  -  sediment  concentration  in  water  column 
r  =  bed  shear  stress 

■*  critical  shear  stress  for  deposition 
Cc  *»  critical  concentration  •=  300  mg /i 


30.  If  the  bed  shear  stress  is  greater  than  the  critical  value  for  par¬ 
ticle  erosion  Tg  ,  material  is  removed  from  the  bed.  The  source  term  is  then 
computed  by  Ariathurai's  (Ariathurai,  MacArthur,  and  Krone  1977)  adaptation  of 
Partheniades '  (1962)  findings: 
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where  P  is  the  erosion  rate  constant,  unless  the  shear  stress  is  also 
greater  than  the  critical  value  for  mass  erosion.  When  this  value  is 
exceeded,  mass  failure  of  a  sediment  layer  occurs  and 
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where 

Tl  -  thickness  of  the  failed  layer 
PL  -  density  of  the  failed  layer 
At  -  time  interval  over  which  failure  occurs 
-  bulk  shear  strength  of  the  layer 

31.  The  cohesive  sediment  bed  consists  of  1  to  10  layers,  each  with  a 
distinct  density  and  erosion  resistance.  The  layers  consolidate  with 
overburden  and  time . 

Bed  shear  stress 

32.  Bed  shear  stresses  are  calculated  from  the  a  ow  speed  according  to 
one  of  four  optional  equations:  the  smooth-wall  log  velocity  profile  or 
Manning  equation  for  flows  alone;  and  a  smooth  bed  or  rippled  bed  equation  for 
combined  currents  and  wind  waves.  Shear  stresses  are  calculated  using  the 
shear  velocity  concept  where 
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where 

=  bed  shear  stress 
u^  -  shear  velocity 

and  the  shear  velocity  is  calculated  by  one  of  four  methods: 
a.  Smooth-wall  log  velocity  profiles 
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which  is  applicable  to  the  lower  15  percent  of  the  boundary 
layer  when 
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where  u  is  the  mean  flow  velocity  (resultant  of  u  and  v 
components) 

b.  The  Manning  shear  stress  equation 
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where  CME  is  a  coefficient  of  1  for  SI  (metric)  units  and 
1.486  for  non- SI  units  of  measurement. 

c.  A  Jonsson-type  equation  for  surface  shear  stress  (plane  beds) 
caused  by  waves  and  currents 
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where 

f  =  shear  stress  coefficient  for  waves 
w 

u  =  maximum  orbital  velocity  of  waves 
om  J 

f  -  shear  stress  coefficient  for  currents 
c 
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d.  A  Bijker-type  equation  for  total  shear  stress  caused  by  waves 
and  current 
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Solution  method 

33.  Equation  B5  is  solved  by  the  finite  element  method  using  Galerkin 
weighted  residuals.  Like  RMA-2V,  which  uses  the  same  general  solution  tech¬ 
nique,  elements  are  quadrilateral  and  may  have  parabolic  sides.  Shape  func¬ 
tions  are  quadratic.  Integration  in  space  is  Gaussian.  Time-stepping  is 
performed  by  a  Crank-Nicholson  approach  with  a  weighting  factor  (0)  of  0.66. 
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A  front- type  solver  similar  to  that  in  RMA-2V  is  used  to  solve  the 
simultaneous  equations. 
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